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INTRODUCTION 

Having  a  narrow  tibia  relative  to  body  mass  has  been  shown  to  be  a  major  predietor  of  stress  fraeture 
risk  and  fragility  (Giladi  et  al,  1987;  Milgrom  et  al,  1989;  Beek  et  al,  1996).  The  reason  for  this 
phenomenon  is  not  fully  understood.  Based  on  studies  of  genetieally  distinet  inbred  mouse  strains,  we 
found  a  reeiproeal  relationship  between  femoral  bone  size  (mid-diaphyseal  eortieal  area)  and  tissue-level 
meehanieal  properties  (bone  quality).  Speeifieally,  we  found  that  slender  bones  (small  width  to  length 
ratio)  were  assoeiated  with  more  damageable  bone  tissue  (Jepsen  et  al,  2001).  We  postulate  that  a 
similar  reciprocal  relationship  between  bone  size  and  bone  material  properties  exists  in  the  human 
skeleton.  The  intriguing  possibility  that  slender  bones,  like  those  we  have  demonstrated  in  animal 
models,  may  be  composed  of  more  damageable  material  than  larger  bones  has  not  been  considered.  To 
test  this  hypothesis,  we  determined  whether  whole  bone  geometry  is  a  predietor  of  tissue  fragility  in  the 
tibiae  from  young  male  donors.  Tissue  damageability  was  assessed  from  biomeehanioal  testing  of 
eompaet  bone  samples  and  eorrelated  with  measures  of  bone  size  and  slenderness.  Speeimens  were 
subjeeted  to  detailed  analyses  of  bone  mierostrueture,  eomposition,  and  mierodamage  eontent.  In  the 
seeond  set  of  experiments,  these  analyses  were  repeated  for  female  donors  to  test  for  gender  differenees 
in  tissue  fragility. 


BODY 

Major  Outcome:  Summary 

All  primary  objeetives  of  this  grant  have  been  sueeessfully  eompleted.  We  found  that  not  all  bones 
are  construeted  in  the  same  manner.  Bones  that  are  more  slender  are  eomprised  of  material  that  is  more 
damageable  eompared  to  bones  that  are  more  robust.  The  inereased  damageability  appears  to  be  a  result 
of  an  inerease  in  ash  eontent.  These  results  were  eonsistent  for  both  males  and  females.  Thus,  young 
adult  males  and  females  with  small  eross-seetional  area  relative  to  bone  length  show  a  eoupled  inerease 
in  tissue-level  mineralization  apparently  in  an  attempt  to  inerease  tissue-stiffness  and  to  eompensate  for 
the  small  bone  geometry  in  order  to  aehieve  an  adequate  whole  bone  stiffness.  The  results  of  this  study 
are  entirely  eonsistent  with  our  investigations  of  the  mouse  skeleton. 

Tissue-level  meehanieal  properties,  ineluding  measures  of  fragility  and  damageability,  were  assessed 
for  eortieal  bone  samples  that  were  maehined  from  the  diaphyses  of  young,  adult  male  and  female  tibiae. 
The  data  indieated  that  the  tibial  diaphyses  of  females  and  males  were  eomposed  of  material  having 
similar  tissue-level  meehanieal  properties,  and  that  the  meehanieal  properties  degraded  with  age  at 
similar  rates  over  the  age  range  of  17-46  years.  Further,  females  and  males  showed  a  similar  relationship 
between  tissue-level  meehanieal  properties  and  eross-seetional  morphology.  These  data  suggested  that 
the  genetie  {Isean,  1984  #491;  Turner,  1990  #492;  Wiren,  2004  #284}  and  environmental  {Gordon, 
1994  #490}  faetors  eontributing  to  sex-speeifie  growth  patterns  affeeted  adult  tibial  eross-seetional  size 
and  shape,  but  did  not  affeet  bone  matrix  eonstruetion,  mineralization,  and  organization  in  a  way  that 
signifioantly  affeeted  tissue-level  meehanieal  properties.  This  analysis  ineluded  meehanieal  properties 
like  stiffness,  whieh  is  important  for  day-to-day  aetivities,  as  well  strength,  duetility,  toughness,  and 
damageability,  whieh  provide  insight  into  the  amount  of  damage  aeeumulated  during  intense  physieal 
aetivity  and  the  fraeture  resistanee  during  an  extreme  load  eondition  sueh  as  a  fall. 

Implications  for  Stress  Fracture  Risk  in  Military  Recruits 

The  results  of  this  study  provide  new  insight  into  why  bone  size  is  a  risk  faetor  for  stress  fraetures. 
Stress  fraetures  are  believed  to  be  a  eonsequenee  of  exeess  damage  aeeumulation  following  intense, 
repetitive  aetivities.  Biologieal  proeesses  that  attempt  to  repair  the  damage  may  further  weaken  the 
tissue  beeause  the  inereased  resorption  results  in  inereased  tissue  porosity  (Sehaffler  and  Burr,  1988). 
However,  the  aetual  eontribution  of  biologieal  repair  proeesses  to  stress  fraeture  risk  remains  unelear 
(Milgrom  et  al,  2004).  Damage,  in  the  form  of  mieroeraeks,  is  the  expeeted  sequelae  of  repetitive 
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loading  following  normal,  daily  activities  (Schaffler  et  al,  1990).  Intense  loading  eonditions,  sueh  as 
those  assoeiated  with  military  training  and  long  distance  running,  are  expeeted  to  further  inerease  in  situ 
damage  aeeumulation  and  degrade  tissue-level  meehanieal  properties  (Mori  et  al,  1993).  Therefore, 
under  extreme  loading  eonditions  (e.g.,  military  training),  variation  in  bone  quality,  speeifieally  tissue 
damageability,  may  be  a  eontributing  faetor  to  the  inereased  risk  of  stress  fraeture  for  individuals  with 
more  slender  bones.  The  eurrent  data  suggested  that  bone  morphology  eould  be  used  as  a  predictor  of 
tissue  fragility  and  stress  fraeture  risk  in  the  absence  of  available  non-invasive  imaging  teehniques  that 
aeeurately  measure  bone  damageability. 

The  results  of  this  study  provided  additional  insight  into  why  female  military  reeruits  show  a  greater 
ineidenee  of  stress  fractures  eompared  to  male  military  reeruits  (Friedl  et  al,  1992;  Beek  et  al  1996).  The 
undersized  morphology  of  female  tibiae  eombined  with  having  similar  tissue-level  meehanieal 
properties  as  males  suggested  that  female  tibiae  may  be  overloaded  during  training  eompared  to  males. 
This  overloading  during  intensive  training  would  be  expeeted  to  lead  to  inereased  tissue-stresses  and 
strains,  and  subsequently  inereased  in  situ  damage  aeeumulation  and  greater  risk  of  developing  a  stress 
fraeture  (Johnson  et  al,  1963;  Burr  et  al,  1990;  Martin  et  al,  1995).  Taken  together,  these  results 
suggested  that  one  faetor  eontributing  to  the  diserepaney  in  stress  fracture  risk  between  females  and 
males  may  be  differenees  in  the  way  bone  size  adapts  to  body  size  after  longitudinal  growth  has  eeased. 
Thus,  fracture  risk  reduetion  would  benefit  from  having  a  better  understanding  of  the  faetors  that 
promote  or  inhibit  bone  adaptation  in  the  young  adult  skeleton  (Lanyon  et  al,  2001;  Milgrom  et  al,  2000; 
Warden  et  al,  2005). 


Major  Outcomes:  Details 
Sample  Population 

We  sueeessfully  aequired  tibiae  from  14  young-adult  female  and  17  young-adult  male  donors.  The 
major  exelusion  eriteria  for  aeeeptance  of  tibia  into  our  study  was  that  the  donors  eould  have  no  medieal 
history  (disease,  pharmaeeutieal  use)  that  would  have  a  major  bone-effeet.  As  shown  in  the  Table  1, 
female  donors  were  shorter  than  males  (163  +  6  em  vs.  178  +  4em;  p<0.0001),  but  were  not  different  in 
body  weight  (74  +  21em  vs.  84  +  25em;  p<0.2). 


Table  1.  Deseription  of  donor  population 


Trait 

Females 

Males 

Number  of  samples 

14 

17 

Average  age 

36.9  ±8.1 

32.9  ±  10.4 

Age  range 

22-46 

17-46 

Weight  (kg) 

73.7±2L1 

83.8  ±25.1 

Height  (em) 

163.5  +  5.9 

177.7  ±4.3 

BMI  (kgW) 

27.4  ±7.2 

26.6  ±8.0 

Tibia  Length  (em) 

33.9  ±2.8 

38.1  ±  1.9 

Monotonie  properties 

Tissue-level  meehanieal  properties  were  assessed  by  loading  four  eortieal  bone  samples  from  each 
tibia  to  failure  in  4-point  bending  at  0.05mm/s  (Fig.  2A)  using  a  servohydraulie  materials  testing  system 
(Instron  model  8872,  Instron  Corp.,  Canton,  MA,  USA).  Speeimens  were  submerged  in  a  PBS  solution 
with  added  ealcium  (Gustafson  et  al,  1996)  and  maintained  at  37°C  throughout  all  tests.  Load  and 
deflection  were  eonverted  to  stress  and  strain  using  the  following  equations  whieh  take  yielding  into 
eonsideration  (Nadai,  1950): 
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(1) 

(2) 


a  =  2[2M  +  (^dM/d(^]/bh^ 

8  =  H/2a  =  1/2  hA[(L  -  a)/(2a^/3  -  a^L  +  lVS)] 

where  a  and  a  are  the  stress  and  strain  at  the  outer  surfaee  of  the  beam,  M  =  applied  moment,  b  = 
speeimen  width,  h  =  speeimen  height,  a  =  14  the  span  between  the  upper  two  load  points  =  9mm,  L  =  14 
the  span  between  the  two  lower  load  points  =  21mm,  and  (j)  =  angle  of  inelination  =  a/p.  The  angle  of 
inelination  was  written  in  terms  of  the  measured  deflection  (A)  by  estimating  the  curvature  (p)  using 
standard  beam  equations.  Mechanical  properties  were  calculated  from  the  stress-strain  curves  and  these 
included  modulus,  strength,  total  energy,  and  post-yield  strain.  Modulus  was  calculated  from  a  linear 
regression  of  the  initial  portion  of  the  stress-strain  curve.  Yield  was  determined  using  the  0.2%  offset 
method.  Post-yield  strain  was  defined  as  the  strain  at  failure  minus  the  strain  at  yield.  All  properties 
were  averaged  over  the  four  samples  tested  for  each  tibia. 

Females  and  males  showed  similar  regressions  for  tissue  modulus,  strength,  post-yield  strain,  and 
energy-to-failure  versus  age  (Figure  1).  No  differences  in  the  average  tissue-level  mechanical  properties 
were  observed  between  females  and  males  (Table  2),  even  when  corrected  for  differences  in  donor  age. 
This  included  measures  of  stiffness  (modulus),  strength,  ductility  (failure  strain,  post-yield  strain),  and 
toughness  (energy-to-failure).  Thus,  females  and  males  show  similar  tissue-level  monotonic  properties. 

Our  data  confirm  that  males  and  females  achieve  similar  adult  tissue-level  mechanical  properties, 
which  has  only  been  assumed  in  prior  analyses  (Ruff,  2000).  Prior  studies  reported  property  values  for 
this  age-range,  but  did  not  directly  compare  data  for  males  and  females  (Currey  and  Butler,  1975; 
Currey,  1979).  The  mechanical  properties  measured  from  the  4-point  bending  tests  in  this  study  were 
consistent  with  the  mechanical  properties  of  cortical  bone  measured  previously  in  tension  (Reilly  and 
Burstein,  1975;  Burstein  et  al,  1976;  McCalden  et  al,  1993),  which  was  expected  given  that  failure  in  our 
bending  tests  occurred  following  fracture  initiated  within  the  tensile  region.  Further,  the  bending 
formula  that  was  used  to  convert  load  and  deflection  to  stress  and  strain,  because  it  accounted  for 
nonlinear  effects  (Nadai,  1950),  provided  a  more  appropriate  estimate  of  bending  strength  compared  to 
standard  beam  theory  (Burstein  et  al,  1972).  Because  the  bone  samples  used  in  this  study  were  machined 
from  within  the  middle  of  the  cortex,  the  tissue  that  was  added  to  the  subperiosteal  or  subendosteal 
surfaces  during  and  after  puberty  was  likely  not  tested  in  our  study.  Thus,  the  current  data  could  not  be 
used  to  determine  if  sex-specific  differences  in  tissue-quality  exist  following  puberty.  This  data  thus 
reflects  the  tissue  properties  established  early  in  life,  plus  the  modifications  to  these  properties 
associated  with  osteonal-remodeling  and  cortical  drift  (Enlow,  1963)  during  the  ensuing  years. 

The  similarity  in  tissue-level  mechanical  properties  for  adult  men  and  women  may  be  explained 
based  on  early  transverse  bone -growth  patterns.  During  pre-pubertal  growth,  the  periosteal  and  endosteal 
surfaces  of  female  and  male  long  bone  diaphyses  follow  nearly  identical  expansion  rates  (Garn,  1970). 
Despite  the  dimorphic  transverse  bone-growth  patterns  arising  during  puberty  (Garn,  1970),  the  long 
bones  of  both  sexes  appear  to  be  constructed  in  a  way  so  that  by  20  years  of  age  bone  size  is  properly 
adapted  to  satisfy  loading  demands  associated  with  body  weight  (Sumner  and  Andriacchi,  1996;  Moro  et 
al,  1996;  Ruff,  2003;  Forwood  et  al,  2004).  Whole  bone  stiffness  and  peak  tissue  strains  are  thus  kept  at 
proper  levels  (Frost,  1987)  during  growth  simply  by  depositing  enough  tissue  on  the  periosteal  and 
endosteal  surfaces  to  keep  pace  with  weight  gain.  Thus,  because  there  is  no  discrepancy  in  the 
relationship  between  bone  size  and  body  size  during  growth,  males  and  females  can  construct  the  tibia 
with  similar  tissue-level  mechanical  properties  without  a  loss  in  mechanical  function. 
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Figure  1.  Females  show  an  age-related  ehange  in  tissue-level  a)  modulus,  b)  strength,  e)  energy-to- 
failure,  and  d)  damageability  that  is  similar  to  males. 


Table  2,  Monotonic  mechanical  properties  for  young-adult  females  and  males. 


MECHANICAL  PROPERTY 

EEMALE 

MALE 

p-value 

Modulus,  GPa 

17.5  ±1.8 

17.1  ±  1.7 

0.5 

Yield  Strain 

0.008  ±  0.0004 

0.008  +  0.0003 

0.7 

Yield  Stress,  MPa 

106.2+10.0 

104.4  +  7.8 

0.6 

Post-yield  Strain 

0.023  +  0.006 

0.025  +  0.005 

0.3 

Eailure  Strain 

0.031  +  0.006 

0.033  +  0.005 

0.3 

Strength,  MPa 

133.8  +  8.4 

130.1+4.2 

0.2 

Energy-to-Eailure,  MPa 

3.1  +  0.7 

3.2 +  0.7 

0.3 

Data  was  age-adjusted  and  shown  as  mean  +  standard  deviation. 


Tissue-Level  Damageability 

Tissue  damageability  was  assessed  using  a  protocol  designed  to  induce  and  accumulate  cracks  in 
cortical  bone  specimens.  The  accumulation  of  damage  leads  to  measurable  degradation  of  mechanical 
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properties  (Lemaitre,  1992).  Therefore,  the  degradation  of  mechanical  properties  can  be  used  as  an  index 
of  matrix  damage.  Four  cortical  bone  samples  from  each  tibia  were  subjected  to  a  fifteen  cycle  damage 
accumulation  protocol  (Figure  2A)  similar  to  that  described  previously  (Jepsen  and  Davy,  1997).  For 
this  protocol,  “diagnostic”  cycles  (1,  3,  5,  7,  9,  11,  13,  and  15)  were  interposed  between  “damage” 
cycles  (2,  4,  6,  8,  10,  12,  and  14).  For  the  diagnostic  cycles,  the  specimens  were  loaded  in  four-point 
bending  at  0.5mm/s  to  50%  of  the  average  displacement  at  yield  (determined  from  the  monotonic  tests), 
held  for  60  seconds,  and  then  unloaded  at  0.5mm/s.  Preliminary  studies  indicated  that  this  load  level 
provided  information  on  tissue-level  mechanical  properties  without  inducing  additional  damage.  For  the 
damage  cycles,  the  specimens  were  loaded  at  0.5mm/s  to  50,  75,  100,  125,  150,  175,  and  200%  of 
displacement  at  yield  respectively,  held  for  60  seconds,  and  then  unloaded  at  0.5mm/sec.  A  5-minute 
recovery  period  followed  each  damage  cycle.  Displacement  at  yield  was  used  as  a  reference  in  the 
damage  cycles  because  this  parameter  showed  little  variation  among  the  test  samples  when  subjected  to 
monotonic  four-point  bending.  The  displacement  at  yield  was  1.0mm  for  the  samples  with  a  height  of 
2.5mm  and  1.07mm  for  the  samples  with  a  height  of  2.2mm. 

A  mechanical  measure  of  the  amount  of  damage  that  accumulated  within  the  test  sample  was 
quantified  from  the  magnitude  of  stiffness  degradation.  For  each  diagnostic  cycle,  stiffness  was 
calculated  from  a  linear  regression  of  the  initial  portion  of  the  load-deformation  curve.  Specimen 
stiffness  decreased  non-uniformly  with  each  cycle  revealing  increasing  amounts  of  damage  induced 
within  each  cycle  and  an  overall  damage  accumulation  by  the  end  of  the  protocol  (Figure  2B).  At  the 
end  of  the  test  sequence,  the  overall  damage  parameter,  D,  was  calculated  by  comparing  the  stiffness  of 
the  first  and  last  diagnostic  tests  such  that: 

D=\-Si5/So,  (3) 

where  S'ls  is  the  stiffness  of  the  last  diagnostic  cycle  and  So  is  the  average  stiffness  of  the  first  two 
diagnostic  cycles  (^i,  5*3)  and  the  first  damage  cycle  (Si). 

No  differences  in  either  total  stiffness  degradation  or  the  total  change  in  relaxation  (Table  3)  were 
observed  between  females  and  males.  Further,  the  changes  in  stiffness  and  relaxation  after  each  damage 
cycle  were  similar  for  females  and  males  (Figure  3).  These  results  indicated  that  the  loading  protocol 
introduced  similar  amounts  of  damage  at  each  load  step  for  females  and  males. 


Table  3.  Comparison  of  total  damage  accumulation  parameters  for  female  and  male  cortical  bone. 


DAMAGE  PARAMETER 

EEMALE 

MALE 

p-value 

Stiffness  Degradation 

0.18  ±0.04 

0.17  ±0.03 

0.6 

Relaxation  Degradation 

1.89  ±0.08 

1.91  ±0.12 

0.6 

Data  was  age-adjusted  and  shown  as  mean  +  standard  deviation. 
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Figure  2,  A)  A  fifteen-cycle  loading  protocol  was  used  to  induce  damage  within  machined  cortical  bone 
specimens  and  to  measure  resultant  stiffness  degradation.  Damage  was  induced  during  cycles  2,  4,  6,  8, 
10,  12,  and  14  by  conducting  relaxation  tests  at  increasing  levels  of  applied  displacement  (expressed 
below  as  a  percentage  of  the  displacement  at  yield).  Diagnostic  cycles  were  interposed  between 
damaging  cycles.  Each  damage  cycle  was  preceded  by  a  diagnostic  cycle  at  50%  of  displacement  at 
yield.  A  5-minute  recovery  period  was  introduced  following  the  damage  cycles  to  relieve  residual 
internal  stresses.  B)  The  change  in  stiffness  calculated  between  sequential  diagnostic  tests  was  plotted 
versus  cycle  number  for  the  damage  tests.  The  dashed  curve  represents  a  specimen  showing  little 
stiffness  degradation  (i.e.,  little  damage  accumulation).  The  solid  line  represents  a  specimen  showing 
large  stiffness  degradation  (i.e.,  more  damage  accumulation). 


A. 


Damage  Parameter:  D  =  1  -  S.i5  /  Sq  D  =  0  (no  damage) 

D  =  1  (fracture) 


Low  D 
High  D 
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Figure  3.  A  plot  of  a)  stiffness  degradation,  Dstiff,  and  b)  relaxation  degradation,  Drelax,  for  each 
cycle  in  the  damage-accumulation  protocol.  Females  and  males  show  similar  curves.  Data  was  not  age- 
corrected. 
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The  Relationship  Between  Morphology  and  Tissue-level  Mechanical  Properties 

Tibia  length  (L)  was  measured  as  the  average  distance  between  the  middle  of  the  talar  trochlear  facet 
and  the  medial  and  lateral  proximal  condyles  (Ruff,  2000)  using  a  large-capacity  slide  caliper  with  an 
accuracy  of  +2. 5mm  (Mantex  Precision,  Haglof  Inc.,  Madison,  MS,  USA).  Tibia  widths  in  the  antero¬ 
posterior  (WidthAp)  and  medial-lateral  (WidthML)  directions  were  measured  using  a  300mm  vernier 
caliper  with  an  accuracy  of  +0.02mm  (Fowler  Company  Inc.,  Newton,  MA,  USA).  Measures  were  taken 
at  10%  intervals  from  30%  to  70%  of  the  total  tibia  length  and  averaged  in  order  to  assess  the  diaphyseal 
morphological  traits  (Miller  and  Purkey,  1980). 

Cross-sectional  morphological  traits  were  quantified  from  3 -mm  thick  mid-diaphyseal  cross-sections 
cut  at  30,  50,  and  70%  of  the  total  tibia  length  using  a  diamond  coated  metallurgical  saw  (Model  660, 
South  Bay  Technology  Inc.,  San  Clemente,  CA,  USA).  Each  cross-section  was  digitally  imaged 
(0.024mm/pixel)  and  analyzed  using  IMAQ  Vision  Builder  (version  6.0,  National  Instruments  Corp., 
Austin,  TX,  USA).  Morphological  traits  included  cortical  area  (CtAr),  the  moments  of  inertia  about  the 
antero-posterior  (Iap)  and  medial-lateral  (Iml)  axes,  the  polar  moment  of  inertia  (J  =  Iap  +  Iml),  and  the 
section  modulus  in  the  AP  (J/WidthAp/2)  and  ML  (J/WidthML/2)  directions.  Moment  of  inertia  and 
section  modulus  were  assessed  because  these  geometric  measures  are  related  to  the  bending  and 
torsional  stiffness  of  intact  tibiae.  A  slenderness  index  (S)  was  calculated  as  the  ratio  of  the  AP  and  ML 
section  modulus  values,  respectively,  to  tibia  length  and  body  weight: 


Sap  - 

(L  *  BW) 

(4) 

(J/(WidthAp/2)) 

Sml  = 

(L  *  BW) 

(5) 

(J/(WidthML/2)) 

where  L  =  tibia  length  (mm)  and  BW  =  body  weight  (kg).  Lor  semantic  reasons,  the  formulae  are  the 
inverse  to  that  used  previously  (Selker  and  Carter,  1989),  such  that  an  increase  in  Sap  or  Sml  indicates  a 
more  slender  (gracile)  bone.  All  morphological  traits  were  averaged  over  the  three  cross-sections  for 
each  tibia. 

None  of  the  morphological  traits  showed  a  significant  correlation  with  age  for  either  males  or 
females  (data  not  shown).  Adult  females  achieved  a  different  bone  size  and  shape  compared  to  males: 
female  tibiae  were  smaller  (J,  CtAr,  Width),  more  slender  (Sap,  Sml,  J/L),  and  showed  a  smaller  J:A 
ratio  compared  to  males  (Table  4).  Lemales  and  males  showed  similar  slopes  for  the  regression  of 
J/WidthAP  versus  body  weight  *  tibial  length  (p<0.17,  ANCOVA).  However,  females  showed  a 
significantly  smaller  intercept  compared  to  males  (p<0.0001,  ANCOVA),  indicating  that  female  tibiae 
were  not  as  well  adapted  to  body  weight  (i.e.,  less  robust)  compared  to  males  (Ligure  4).  Similar  results 
were  observed  when  using  J/WidthML- 

The  tissue-level  mechanical  properties  were  regressed  against  bone  morphological  traits  and  no 
significant  differences  in  slopes  were  found  for  females  and  males.  A  correlation  analysis  (Table  5), 
which  was  conducted  using  the  combined  female  and  male  data  sets,  showed  that  tissue-level  modulus 
decreased  with  increasing  bone  section  modulus  (J/WidthAp,  JAVidthMp),  post-yield  strain  increased  with 
increasing  bone  width  (WidthAp),  and  tissue-damage  (Dstiff)  increased  with  increasing  tibial  slenderness 
(Sap,  Sml).  The  significant  regressions  are  shown  in  Ligure  5.  The  amount  of  overlap  between  the 
female  and  male  datasets  depended  on  whether  the  bone  morphology  traits  were  adjusted  for  body  size. 
Lor  unadjusted  traits  like  width  and  J/width,  the  female  data  simply  extended  the  relationships  observed 
for  the  males  (Ligures  5A,  5B).  Lor  the  regression  involving  bone  slenderness,  which  is  adjusted  for 
body  weight  and  tibial  length,  the  female  data  overlapped  substantially  with  the  male  data  (Ligure  5C). 

Although  each  of  the  tissue-level  mechanical  properties  that  we  measured  showed  a  small  coefficient 
of  variation  for  both  sexes,  approximately  20%  of  this  variation  was  explained  by  the  size  of  the  tibia. 
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Female  tibiae  extended  the  eorrelations  that  were  observed  previously  for  males  (Tommasini  et  al, 
2005),  and  indieated  that  individuals  with  smaller  eross-seetional  size  (i.e.,  more  slender  bones)  showed 
a  larger  tissue-level  stiffness.  This  relationship  may  be  indieative  of  an  inherent  adaptive  response  of 
bone  to  modify  tissue-level  meehanieal  properties  in  order  to  eompensate  for  the  smaller  size  (Ferretti  et 
al,  1993).  This  eoupling  suggested  that  osteoblasts  and  osteoelasts  are  eapable  of  eoordinately  adapting 
bone  morphology  and  tissue-level  meehanieal  properties  so  that  the  eombination  of  these  traits  satisfies 
meehanieal  demands.  A  more  dramatic  coordinated  relationship  among  morphological  traits,  tissue- 
quality,  and  mechanical  function  has  been  observed  previously  across  a  large  range  of  bones  from  many 
species  (Currey,  1979)  and  for  inbred  mice  (Jepsen  et  al,  2001;  Tommasini  et  al,  2005;  Price  et  al, 
2005),  bats  (Swartz  et  al,  1992),  gulls  (Carrier  and  Leon,  1990),  and  polar  bears  (Brear  et  al,  1990). 
Although  this  coupling  appears  advantageous  for  ensuring  that  an  adequate  whole  bone  stiffness  is 
achieved  for  day-to-day  activity,  the  disadvantage  is  that  the  tissue-quality  factors  that  tend  to  make 
bone  stiff  also  tend  to  make  bone  less  ductile,  less  tough,  and  more  damageable  (Currey,  1984).  These 
latter  tissue-quality  factors  are  important  during  extreme  loading  conditions  (Turner,  2002). 

The  relationship  between  bone  size  (section  modulus)  and  body  size  (body  weight  *  tibia  length) 
was  established  previously  by  Ruff  (Ruff,  2000)  and  was  based  on  the  resistance  of  bone  to  bending 
loads.  The  divergence  between  males  and  females  suggested  that  the  adaptive  factors  that  match  bone 
size  to  body  size  during  growth  (Sumner  and  Andriacchi,  1996;  Moro  et  al,  1996;  Ruff,  2003;  Forwood 
et  al,  2004)  no  longer  work  the  same  way  when  men  and  women  reach  adulthood.  This  latent  difference 
in  bone  adaptation  may  be  due  to  differences  in  post-pubescent  endocrine  factors  that  regulate  the 
relative  amounts  of  bone  added  to  the  periosteum  versus  the  endosteum  (Duan  et  al,  2003).  Our  data, 
and  those  of  others  (Geusens  et  al,  1991;  Looker  et  al,  2001;  Duan  et  al,  2003;  Nieves  et  al,  2005), 
indicated  that  the  relationship  between  bone  size  and  body  size  diverges  for  females  and  males  after  20 
years  of  age  (Figure  4),  such  that  female  long  bones  become  progressively  undersized  relative  to  body 
size.  Thus,  adult  females  would  be  expected  to  accumulate  more  damage  under  intense  loading 
compared  to  males  and  this  may  be  a  contributing  factor  to  the  larger  stress  fracture  incidence  for  female 
military  recruits  (Friedl  et  al,  1992;  Beck  et  al,  1996). 


Table  4,  Diaphyseal  cross-sectional  morphology  for  female  and  male  tibiae. 


MORPHOLOGICAL  TRAIT 

FEMALE 

MALE 

p-value 

TIBIA  LENGTH,  mm 

33.9  +  2.8 

38.1  +  1.9 

0.0001 

WIDTHap,  mm 

26.5  +  2.7 

31.2  +  2.5 

0.0001 

WIDTHml,  mm 

20.8  +  2.4 

24.3  +  2.3 

0.0001 

CtAr,  mm^ 

249  ±  40 

356  ±55 

0.0001 

J,  mm 

24776  +  7917 

51640  ±  15886 

0.0001 

J  /  rAP,  m^ 

1836  +  473 

3279  ±819 

0.0001 

J  /  pml, 

2352  ±613 

4188  ±907 

0.0001 

J  /  L,  m^ 

722 ± 194 

1352  ±405 

0.0001 

J  /  A,  m^ 

97.1  ±  18.5 

142.7  ±24.9 

0.0001 

Sap,  1/mm  /kg 

14.1+4.1 

9.9  ±2.0 

0.002 

Sml,  1/mm^/kg 

11.1  +  3.6 

7.7  ±  1.5 

0.004 

Data  was  calculated  over  30-70%  of  the  total  tibial  length  and  shown  as  mean  +  standard  deviation. 
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Table  5.  Correlation  matrix  between  morphologieal  traits  and  tissue-level  meehanieal  properties. 


MODULUS 

STRENGTH 

PY  STRAIN 

ENERGY 

DAMAGE 

WIDTHap 

-0.13 

-0.15 

0.35 

0.34 

-0.12 

(0.5) 

(0.4) 

(0.05) 

(0.06) 

(0.5) 

WIDTHml 

-0.13 

-0.17 

0.29 

0.28 

-0.03 

(0.5) 

(0.4) 

(0.12) 

(0.13) 

(0.8) 

CtAr 

-0.25 

-0.18 

0.31 

0.30 

-0.24 

(0.2) 

(0.3) 

(0.09) 

(0.1) 

(0.2) 

J 

-0.33 

-0.25 

0.27 

0.24 

-0.24 

(0.07) 

(0.2) 

(0.14) 

(0.2) 

(0.2) 

J/L 

-0.35 

-0.26 

0.26 

0.22 

-0.26 

(0.05) 

(0.2) 

(0.2) 

(0.2) 

(0.2) 

J  /  WidthAP 

-0.36 

-0.25 

0.22 

0.18 

-0.27 

(0.04) 

(0.2) 

(0.2) 

(0.3) 

(0.1) 

J  /  WidthML 

-0.36 

-0.26 

0.24 

0.20 

-0.30 

(0.05) 

(0.2) 

(0.2) 

(0.3) 

(0.1) 

Sap 

0.25 

-0.05 

-0.10 

-0.09 

0.43 

(0.2) 

(0.8) 

(0.6) 

(0.7) 

(0.02) 

C 

0.27 

-0.04 

-0.13 

-0.11 

0.44 

OML 

(0.1) 

(0.8) 

(0.56) 

(0.6) 

(0.01) 

Combined  male  and  female  data  sets. 

First  row  =  Pearson  correlation  coefficient.  Second  row  =  p-value. 
Damage  =  stiffness  degradation. 


Figure  4.  Variation  in  the  section  modulus  (Polar  moment  of  inertia,  J  /  WidthAp)  versus  body  size 
(body  weight  *  tibial  length)  for  adult  female  and  male  tibiae. 


BODY  WEIGHT  X  BONE  LENGTH  (kg  mm) 
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Figure  5.  Variation  in  mechanical  properties  as  a  function  of  diaphyseal  cross-sectional  morphology. 
Solid  lines  represent  the  regression  for  the  combined  female  and  male  datasets.  Dashed  lines  represent 
the  regressions  for  males  and  females  separately. 


0.0  2.5  5.0  7.5  10.0  12.5  15.0  17.5  20.0  22.5 


SLENDERNESSap  (t/mm^/kg) 
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The  Relationship  Among  Morphological,  Tissue-level  Mechanical  Properties,  Tissue  Composition,  and 
Matrix  Architecture 

The  density,  ash  content,  and  water  content  were  determined  for  each  sample  retrieved  from  the 
monotonic  tests.  Specimen  volume  was  determined  using  Archimedes’  principal.  Submerged  weight, 
hydrated  weight,  dry  weight,  and  ash  weights  were  measured.  To  test  for  variation  in  matrix 
organization  in  males  and  females,  bone  microstructure  was  assessed  for  each  sample  retrieved  from  the 
damageability  tests.  For  each  specimen,  digital  images  of  three  transverse  sections,  lOOpm  in  thickness, 
were  taken  at  lOX,  stitched  together,  and  traced  using  a  Tablet  monitor  (WACOM).  Parameters 
measured  include  porosity  and  the  area  fractions  of  osteonal,  interstitial  (remodeled),  and 
circumferential  lamellar  (unremodeled)  tissues.  Both  vascular  canals  and  resorption  spaces  were  counted 
as  pores.  Osteonal  tissue  was  defined  as  a  lamellar  region  with  a  Haversian  canal  completely  surrounded 
by  a  cement  line.  Data  from  individual  test  samples  were  averaged  for  each  donor. 

To  determine  if  males  and  females  show  a  similar  relationship  between  tibia  cross-sectional 
morphology,  bone  microstructure,  and  tissue-level  mechanical  properties,  linear  regressions  were 
performed  while  taking  age  into  consideration  and  slopes  and  intercepts  between  male  and  female 
regressions  were  compared  (ANCOVA). 

Since  no  differences  were  found  between  male  and  female  regressions,  combined  male  and  female 
datasets  were  used  to  compare  tissue-level  mechanical  properties,  cross-sectional  morphology, 
composition,  and  bone  microstructure.  The  only  significant  correlations  between  compositional  or 
architectural  traits  and  tissue-level  mechanical  properties  were  observed  between  ash  content  and  tissue 
modulus  (p  <  0.003),  strength  (p  <  0.02),  and  post-yield  strain  (p  <0.01;  Figure  6A).  Measures  of  bone 
size  (total  cross-sectional  area,  cortical  area,  ML  and  AP  width,  cross-sectional  polar  moment  of  inertia) 
were  negatively  correlated  with  ash  content  (p  <  0.02).  Bone  slenderness  (polar  moment  of  inertia 
normalized  for  width  and  polar  moment  of  inertia  normalized  for  tibia  length)  was  also  negatively 
correlated  with  ash  content  (p  <  0.02;  Figure  6B).  Bone  size  and  slenderness  were  positively  correlated 
with  the  area  fraction  of  remodeled  tissue  (p  <  0.05).  Thus,  wider  (more  robust)  tibiae  had  lower  mineral 
content  and  increased  amount  of  remodeled  tissue  and  narrow  (more  slender)  tibiae  had  increased 
mineral  content  and  decreased  amount  of  remodeled  tissue. 

The  current  study  confirmed  that,  similar  to  inbred  mice,  the  tibia  diaphyses  of  young,  adult  females 
and  males  show  a  reciprocal  relationship  between  bone  morphology  and  tissue  quality,  as  hypothesized. 
Tissue  from  more  slender  bones  (narrow  relative  to  tibia  length)  was  composed  of  material  having 
higher  mineral  content  and,  consequently,  increased  tissue  stiffness  and  strength,  but  reduced  ductility. 
Unlike  the  mouse  model,  this  relationship  was  subtle  and  only  observed  when  combining  male  and 
female  datasets  to  expand  the  range  of  data.  The  current  study  revealed  a  general  biological  phenomenon 
that  has  only  been  observed  previously  by  looking  across  a  large  range  of  bones  from  many  species. 
Here,  bone  from  the  same  species  with  similar  function  showed  that  more  slender  morphology  might  be 
compensated  by  increased  ash  content.  Although  this  coupling  appears  advantageous  for  ensuring  that 
an  adequate  whole  bone  stiffness  is  achieved  for  day-to-day  activity,  the  disadvantage  is  that  the  tissue- 
quality  factors  that  tend  to  make  bone  stiff  also  tend  to  make  bone  less  ductile,  less  tough,  and  more 
damageable. 

This  is  an  intriguing  finding  because,  despite  dimorphic  growth  patterns,  males  and  females  have 
similar  tissue-level  mechanical  properties.  Thus,  both  sexes  may  share  a  common  underlying  biological 
(adaptive)  mechanism.  This  data  suggest  that  osteoblasts  and  osteoclasts  are  capable  of  not  only 
adapting  bone  morphology,  but  also  modulating  tissue-level  mechanical  properties  so  that  the 
combination  of  traits  satisfies  mechanical  demands  (Jepsen  et  al,  2001;  Currey,  1979).  Therefore,  the 
phenomenon  that  slender  bones  are  constructed  differently  with  material  level  variation  that  ultimately 
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leads  to  more  damageable  material  than  larger  bones  may  help  explain  why  individuals  with  more 
slender  bones  are  at  inereased  risk  of  stress  fractures  early  in  life  and  fragility  fractures  later  in  life. 


Figure  6.  Correlation  between  ash  content  and  (A)  post-yield  strain,  PYs  and  (B)  polar  moment  of 
inertia,  J,  normalized  for  tibia  length,  L. 


Ultrasound  Studies 

In  year  three,  we  conducted  a  series  of  experiments  that  examined  whether  ultrasound  has  the 
sensitivity  to  detect  the  presence  of  damage  in  bone.  Most  studies  have  found  that  ultrasound  is 
insensitive  to  damage.  We  found  that  damage  affects  the  viscoelastic  properties  of  bone  to  a  greater 
degree  than  the  stiffness  properties  and  we  focused  our  attention  on  the  attenuation  of  the  ultrasound 
signal  as  this  may  be  more  sensitive  to  the  presence  of  damage.  A  student  for  the  City  University  of  New 
York  worked  on  this  project  and  found  that,  like  prior  work,  the  velocity  of  the  ultrasound  signal  was 
insensitive  to  the  presence  of  damage.  He  examined  the  attenuation  (a  measure  related  to  the  viscoelastic 
properties  of  bone)  following  damage  accumulation,  however,  could  not  make  a  reasonable  conclusion 
because  of  difficulties  in  matching  the  size  of  the  sample  to  the  size  of  the  ultrasound  transducer  (i.e.,  we 
would  have  to  use  samples  machined  from  bovine  bone  to  test  this  concept).  Although  viscous 
properties  may  be  more  sensitive  to  the  presence  of  damage,  we  concluded  from  this  study  that  it  was 
not  practical  to  reliably  measure  this  parameter  using  machined  bone  samples  from  human  tibiae.  Thus, 
it  is  uncertain  whether  ultrasound  would  be  a  valuable  measure  to  test  for  the  presence  of  in  situ  damage 
accumulation  engendered  during  basic  training. 
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KEY  RESEARCH  ACCOMPLISHMENTS 

The  primary  outcome  of  the  grant  is  that  the  bone  tissue  for  individuals  with  more  slender  bones  are 
constructed  in  a  fundamentally  different  way  compared  to  individuals  with  more  robust  (wider)  bones. 
With  increasing  slenderness,  the  material  is  more  stiff,  less  ductile,  and  more  damageable  and  this 
appears  to  be  a  consequence  of  higher  matrix  mineralization.  These  results  support  are  central 
hypothesis  and  may  help  explain  why  individuals  with  smaller  tibiae  are  at  higher  risk  of  stress 
fractures. 
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Funding 

The  data  generated  by  this  grant  provided  evidence  that  the  mouse  represents  an  important  model  for 
understanding  the  genetic  variability  in  the  human  skeleton.  This  data  helped  secure  an  ROl  grant  from 
the  NIH  (AR44927)  titled,  "Genetic  Determination  of  Skeletal  Fragility". 

The  data  generated  by  this  grant  was  also  used  in  the  submission  of  two  NIH  grants  titled,  “Genetic 
regulation  of  bone  growth  and  development”  and  “Age-changes  in  bone  morphology  as  a  predictor  of 
fracture  risk”.  The  goal  of  the  first  grant  is  to  identify  chromosomes  harboring  quantitative  trait  loci 
influencing  bone  growth  patterns.  The  goal  of  the  second  grant  is  to  determine  if  25  year  changes  in 
moment  of  inertia  predict  fracture  incidence.  This  latter  grant  utilizes  information  from  the  Framingham 
Heart  Study  and  was  written  in  collaboration  with  the  Boston  University  School  of  Public  Health  and  the 
Hebrew  Rehabilitation  Center  for  the  Aged. 


CONCLUSIONS 

The  results  of  this  grant  have  revealed  a  novel  relationship  between  bone  morphology  and  tissue 
mechanical  properties.  The  investigations  of  the  mouse  skeleton  revealed  that  genetic  variation  in  bone 
morphology  strongly  influence  tissue  mechanical  properties  through  variations  in  matrix  composition. 
The  data  from  this  DOD  grant  indicated  that  a  similar  relationship  also  exists  in  the  human  skeleton. 
Thus,  individuals  who  have  smaller  (more  narrow)  tibia  for  their  body  size  appear  to  compensate  for  the 
smaller  geometry  through  variation  in  tissue-level  mechanical  properties.  One  of  the  side  effects  of  this 
compensation  is  altered  damageability  which  may  be  revealed  under  extreme  physical  activity  such  as 
that  experienced  during  military  training. 
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Relationship  Between  Bone  Morphology  and  Bone  Quality  in  Male 
Tibias:  Implications  for  Stress  Fracture  Risk 


Steven  M  Tommasini,^  Philip  Nasser,^  Mitchell  B  Schaffler,^  and  Karl  J  Jepsen^ 


ABSTRACT:  Biomechanical  properties  were  assessed  from  the  tihias  of  17  adult  males  17-46  years  of  age. 
Tissue-level  mechanical  properties  varied  with  hone  size.  Narrower  tihias  were  comprised  of  tissue  that  was 
more  brittle  and  more  prone  to  accumulating  damage  compared  with  tissue  from  wider  tihias. 

Introduction:  A  better  understanding  of  the  factors  contributing  to  stress  fractures  is  needed  to  identify  new 
prevention  strategies  that  will  reduce  fracture  incidence.  Having  a  narrow  (i.e.,  more  slender)  tibia  relative  to 
body  mass  has  been  shown  to  be  a  major  predictor  of  stress  fracture  risk  and  fragility  in  male  military  recruits 
and  male  athletes.  The  intriguing  possibility  that  slender  bones,  like  those  shown  in  animal  models,  may  be 
composed  of  more  damageable  material  has  not  been  considered  in  the  human  skeleton. 

Materials  and  Methods:  Polar  moment  of  inertia,  section  modulus,  and  antero-posterior  (AP)  and  medial- 
lateral  (ML)  widths  were  determined  for  tibial  diaphyses  from  17  male  donors  17^6  years  of  age.  A  slen¬ 
derness  index  was  defined  as  the  inverse  ratio  of  the  section  modulus  to  tibia  length  and  body  weight.  Eight 
prismatic  cortical  bone  samples  were  generated  from  each  tibia,  and  tissue-level  mechanical  properties  in¬ 
cluding  modulus,  strength,  total  energy,  postyield  strain,  and  tissue  damageability  were  measured  by  four- 
point  bending  from  monotonic  (n  =  4/tibia)  and  damage  accumulation  {n  =  4/tibia)  test  methods.  Partial 
correlation  coefficients  were  determined  between  each  geometrical  parameter  and  each  tissue-level  mechani¬ 
cal  property  while  taking  age  into  consideration. 

Results:  Significant  correlations  were  observed  between  tibial  morphology  and  the  mechanical  properties  that 
characterized  tissue  brittleness  and  damageability.  Positive  correlations  were  observed  between  measures  of 
bone  size  (AP  width)  and  measures  of  tissue  ductility  (postyield  strain,  total  energy),  and  negative  correlations 
were  observed  between  bone  size  (moment  of  inertia,  section  modulus)  and  tissue  modulus. 

Conclusions:  The  correlation  analysis  suggested  that  bone  morphology  could  be  used  as  a  predictor  of  tissue 
fragility  and  stress  fracture  risk.  The  average  mechanical  properties  of  cortical  tissue  varied  as  a  function  of  the 
overall  size  of  the  bone.  Therefore,  under  extreme  loading  conditions  (e.g.,  military  training),  variation  in  bone 
quality  parameters  related  to  damageability  may  be  a  contributing  factor  to  the  increased  risk  of  stress  fracture 
for  individuals  with  more  slender  bones. 
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INTRODUCTION 

Stress  fractures  are  overuse  injuries  of  bone  that  are 
common  among  elite  runners  and  military  recruits.^^^^^ 
Before  injury,  affected  bones  are  typically  normal  with  no 
acute  injury.  Morbidity  from  stress  fractures  ranges  from 
minor  pain  to  serious  lifetime  disability  for  the  individual.*^"*^ 
Stress  fractures  have  been  reported  in  the  ribs,  hip,  spine, 
and  metatarsals, but  vigorous  weight-bearing  activities, 
such  as  running  and  jogging,  commonly  lead  to  stress  frac¬ 
tures  of  the  lower  extremities,  especially  the  tibia.^^^  During 
basic  training,  1-5%  of  U.S.  male  military  recruits  sustain  a 
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stress  fracture.^^^  However,  this  incidence  is  two  to  five 
times  higher  in  female  recruits. Stress  fractures  lead  to 
loss  of  manpower,  valuable  loss  of  training  time,  expense  of 
medical  care,  and  discharge  of  affected  soldiers. A  better 
understanding  of  the  factors  contributing  to  stress  fractures 
is  needed  to  identify  new  prevention  strategies  that  will 
reduce  fracture  incidence. 

A  number  of  risk  factors  for  stress  fracture  have  been 
identified  including  physical  fitness,  external  hip  rotation, 
body  height  and  weight,  age,  race,  gender,  muscle  mass, 
motivation,  footwear,  smoking,  and  family  history  of  osteo¬ 
porosis. One  of  the  best  predictors  of  stress  fracture 
risk  is  bone  geometry.  Specifically,  having  a  narrow  (i.e., 
more  slender)  tibia  relative  to  body  mass  has  been  shown  to 
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1372 


RELATIONSHIP  BETWEEN  BONE  MORPHOLOGY  AND  QUALITY 


1373 


be  a  major  predictor  of  stress  fracture  risk  and  fragility  in 
male  military  recruits^^’^’^^^  and  male  athletes. A  stress 
fracture  is  thought  to  be  a  consequence  of  transiently  re¬ 
duced  tissue  strength  arising  from  increased  resorptive  ac¬ 
tivity  (i.e.,  increased  porosity)  that  acts  to  repair  damage 
induced  by  vigorous  physical  activity.^^^^  Thus,  stress  frac¬ 
tures  may  be  pronounced  in  individuals  with  more  slender 
bones  because  smaller  bone  size  is  thought  to  lead  to  higher 
tissue-level  stresses  and  thus  increased  damage  accumula- 
tion.*-^'^^  However,  this  postulate  is  based  on  the  assumption 
that  all  bones  are  constructed  in  equivalent  manners,  and 
the  contribution  of  variable  tissue -level  mechanical  proper¬ 
ties  to  stress  fracture  incidence  has  not  been  explored. 

An  examination  of  inbred  mouse  strains  may  help  ex¬ 
plain  why  bone  size  is  a  risk  factor  for  stress  fractures  in  the 
human  skeleton.  A  comparison  of  adult  A/J  and  C57BL/6J 
inbred  mouse  strains  revealed  that  the  bone  slenderness 
was  inversely  related  to  mineral  content  (as  measured  by 
ash  content)  and,  by  correlation,  tissue  modulus  and 
strength.^^"*^  Mineral  content  has  been  shown  to  be  posi¬ 
tively  correlated  with  tissue  stiffness  and  strength.^^^^  These 
results  suggested  that  bone  morphology  and  mineral  con¬ 
tent  were  coordinately  regulated  so  whole  bone  stiffness 
appropriately  matched  the  mechanical  demands  imposed 
by  weight  bearing.  However,  the  downside  of  regulating 
mineral  content  to  match  bone  size  was  that  mineral  con¬ 
tent  was  also  negatively  correlated  with  tissue  ductil- 
ity  (14.15)  postulate  that  a  similar  reciprocal  relationship 
between  bone  size  and  bone  quality  exists  in  the  human 
skeleton.  The  intriguing  possibility  that  slender  bones,  like 
those  shown  in  animal  models,  may  be  composed  of  more 
damageable  material  has  not  yet  been  considered  in  the 
human  skeleton. 

The  goal  of  this  study  was  to  determine  whether  tissue- 
level  mechanical  properties  vary  with  bone  size  in  the  human 
skeleton.  This  was  tested  by  assessing  the  biomechanical 
properties  of  tibias  from  young  adult  males.  Understand¬ 
ing  why  bone  morphology  is  a  risk  factor  for  stress  frac¬ 
tures  should  lead  to  better  identification  of  those  at  risk 
and,  ultimately,  to  early  diagnosis,  treatment,  and  modifi¬ 
cation  of  training  regimens. 

MATERIALS  AND  METHODS 

Sample  population 

Tibias  of  17  male  donors  (15  white,  1  Hispanic,  1  black) 
32.9  ±  10.4  years  of  age  (range,  17^6  years)  were  acquired 
from  the  Musculoskeletal  Transplant  Foundation  (Edison, 
NJ,  USA).  Donor  body  weight  and  height  were  obtained 
from  the  source.  Only  donors  with  no  known  skeletal  pa¬ 
thology  were  included  in  the  study.  The  tibias  were  freshly 
harvested,  wrapped  in  wet  gauze,  and  stored  in  plastic  bags 
at  ^0°C. 

Whole  bone  morphology 

Tibia  length  (L)  was  measured  as  the  average  distance 
between  the  distal  articular  center  (the  middle  of  the  talar 
trochlear  facet)  and  the  two  proximal  articular  centers  (me¬ 
dial  and  lateral  condyles)*^^®^  using  a  large-capacity  slide 
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FIG.  1.  Schematic  of  how  whole  tibias  were  sectioned  to  pro¬ 
duce  the  3-mm-thick  sections  used  for  cross-sectional  morphology 
and  cortical  bone  samples  for  biomechanical  testing  from  three 
diaphyseal  cylindrical  sections  (monotonic,  n  =  4/tibia;  damage 
accumulation,  n  =  4/tibia). 


caliper  with  an  accuracy  of  ±2.54  mm  (Mantex  Precision; 
Haglof,  Madison,  MS,  USA).  Tibia  width  was  measured  in 
the  antero-posterior  (AP)  and  medial-lateral  (ML)  direc¬ 
tions  at  10%  intervals  from  30-70%  of  the  total  tibia  length 
using  a  300-mm  vernier  caliper  with  an  accuracy  of  ±0.02 
mm  (Fowler  Company,  Newton,  MA,  USA). 

Cross-sectional  morphology  was  determined  from  3-mm- 
thick  middiaphyseal  cross-sections  cut  at  30%,  50%,  and 
70%  of  the  total  tibia  length  (Fig.  1)  using  a  diamond 
coated  metallurgical  saw  (Model  660;  South  Bay  Technol¬ 
ogy,  San  Clemente,  CA,  USA).  A  calibrated  image  of  each 
cross-section  was  obtained  using  a  digital  camera  at  a  0.024 
mm/pixel  resolution.  Image  analysis  software  (IMAQ  Vi¬ 
sion  Builder  6.0;  National  Instruments,  Austin,  TX,  USA) 
was  used  to  threshold  each  image  and  quantify  cortical  area 
(CtAr),  the  moments  of  inertia  about  the  AP  (/^p)  and  ML 
(7^^l)  axes,  the  polar  moment  of  inertia  (J  =  /^p  + 
and  the  section  modulus  in  the  AP  (J/APwidth/2)  and  ML 
(J/MLwidth/2)  directions.  Moment  of  inertia  and  section 
modulus  were  assessed  because  these  geometric  measures 
are  related  to  the  bending  and  torsional  stiffness  of  intact 
tibias.  A  slenderness  index  (S)  was  calculated  in  the  AP  and 
ML  directions  as  the  ratio  of  the  AP  and  ML  section  modu¬ 
lus  values,  respectively,  to  tibia  length  and  body  weight^^^^: 

S  =  l/{[J/(width/2)]/(L  X  BW))  (1) 

where  L  =  tibia  length  (mm)  and  BW  =  body  weight  (kg). 
The  section  modulus  has  been  shown  to  scale  linearly  with 
body  mass.^^^^  The  inverse  ratio  was  used  so  that  a  tibia 
with  a  large  slenderness  value  is  one  that  is  thinner  or  grac- 
ile  for  the  weight  and  height  of  an  individual.  A  small  slen¬ 
derness  value  reflects  a  stocky  or  robust  tibia.  All  morpho¬ 
logical  traits  were  averaged  over  the  three  cross-sections  for 
each  tibia. 

Bone  sample  generation 

Cortical  bone  samples  were  prepared  from  the  diaphysis 
of  each  tibia  for  biomechanical  testing  (Fig.  1).  The  three 
diaphyseal  cylindrical  sections  were  rough-cut  into  antero- 
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FIG.  2.  (A)  Schematic  of  the  four-point  bending  device  shows 

how  cortical  bone  samples  were  tested.  The  span  between  the 
upper  contact  points  is  2a  =  18  mm,  and  the  span  between  the 
lower  contact  points  is  2L  =  42  mm.  (B)  Typical  stress-strain 
curve  from  the  four-point  bending  monotonic  tests  shows  how 
modulus,  strength,  postyield  strain,  and  total  energy  were  calcu¬ 
lated. 


lateral,  antero-medial,  and  posterior  regions.  From  each  of 
these  regions,  one  to  three  prismatic  beams  were  cut  using 
a  diamond-coated  metallurgical  saw  (Isomet;  Buehler, 
Lake  Bluff,  IL,  USA).  The  beams  were  machined  to  regular 
test  samples  using  an  automated  CNC  milling  machine  un¬ 
der  constant  irrigation  (Modela  MDX-20;  Roland  DGA, 
Irvine,  CA,  USA).  Sample  width  (circumferential  direc¬ 
tion)  was  machined  to  5  mm  and  length  (longitudinal  di¬ 
rection)  was  machined  to  55  mm  for  all  samples.  Sample 
height  (radial  direction)  was  2.5  mm  except  for  four  tibias 
with  thin  cortices,  which  were  machined  to  2.2  mm.  A  total 
of  eight  samples  were  generated  from  each  tibia  and  ran¬ 
domly  distributed  to  monotonic  (n  =  4)  and  damage  accu¬ 
mulation  (n  =  4)  test  groups.  All  samples  were  stored  at 
^0°C  in  gauze  saturated  with  PBS  with  added  calcium^^^^ 
and  placed  individually  in  airtight  bags. 

Monotonic  failure  properties 

Tissue-level  mechanical  properties  were  assessed  by 
loading  four  cortical  bone  samples  from  each  tibia  to  failure 
in  four-point  bending  at  0.05  mm/s  (Fig.  2A)  using  a  servo- 
hydraulic  materials  testing  system  (Instron  model  8872;  In- 
stron.  Canton,  MA,  USA).  Specimens  were  submerged  in  a 
PBS  solution  with  added  calcium^^^^  and  maintained  at 
37°C  throughout  all  tests.  Load  and  deflection  were  con¬ 
verted  to  stress  and  strain  using  the  following  equations, 
which  take  yielding  into  consideration*^^®^: 


a  =  2[2M  -t  <|)(dM/d<|))]/bh^  (2) 

E  =  h(|)/2a  =  y2hA[(L  -  a)/(2a^/3  -  a^L  -t  L^/3)]  (3) 

where  a  and  e  are  the  stress  and  strain  at  the  outer  surface 
of  the  beam,  M  =  applied  moment,  b  =  specimen  width,  h 
=  specimen  height,  a  =  1/2  the  span  between  the  upper  two 
load  points  =  9  mm,  L  =  Vi  the  span  between  the  two 
lower  load  points  =  21  mm,  (f*  =  angle  of  inclination  =  a/p, 
and  d/d4>  is  the  derivative  with  respect  to  4>.  The  angle  of 
inclination  was  written  in  terms  of  the  measured  deflec¬ 
tion  (A)  by  estimating  the  curvature  (p)  using  standard 
beam  equations.  Mechanical  properties  were  calculated 
from  the  stress-strain  curves,  and  these  included  modulus, 
strength,  total  energy,  and  postyield  strain  (Fig.  2B).  Modu¬ 
lus  was  calculated  from  a  linear  regression  of  the  initial 
portion  of  the  stress-strain  curve.  Yield  was  determined 
using  the  0.2%  offset  method.  Postyield  strain  was  defined 
as  the  strain  at  failure  minus  the  strain  at  yield.  All  prop¬ 
erties  were  averaged  over  the  four  samples  tested  for  each 
tibia. 

Damage  accumulation  tests 

Tissue  damageability  was  assessed  using  a  protocol  de¬ 
signed  to  induce  and  accumulate  cracks  in  cortical  bone 
specimens.  The  accumulation  of  damage  leads  to  measur¬ 
able  degradation  of  mechanical  properties.^^°^  Therefore, 
the  degradation  of  mechanical  properties  can  be  used  as  an 
index  of  matrix  damage.  Four  cortical  bone  samples  from 
each  tibia  were  subjected  to  a  fifteen  cycle  damage  accu¬ 
mulation  protocol  (Fig.  3A)  similar  to  that  described  pre- 
viously.*^^^^  For  this  protocol,  “diagnostic”  cycles  (1,  3,  5,  7, 
9, 11, 13,  and  15)  were  interposed  between  “damage”  cycles 
(2,  4,  6,  8,  10,  12,  and  14).  For  the  diagnostic  cycles,  the 
specimens  were  loaded  in  four-point  bending  at  0.5  mm/s  to 
50%  of  the  average  displacement  at  yield  (determined  from 
the  mono  tonic  tests),  held  for  60  s,  and  unloaded  at  0.5 
mm/s.  Preliminary  studies  indicated  that  this  load  level  pro¬ 
vided  information  on  tissue-level  mechanical  properties 
without  inducing  additional  damage.  For  the  damage 
cycles,  the  specimens  were  loaded  at  0.5  mm/s  to  50%,  75%, 
100%,  125%,  150%,  175%,  and  200%  of  displacement  at 
yield,  respectively,  held  for  60  s,  and  unloaded  at  0.5  mm/s. 
A  5-minute  recovery  period  followed  each  damage  cycle. 
Displacement  at  yield  was  used  as  a  reference  in  the  dam¬ 
age  cycles  because  this  parameter  showed  little  variation 
among  the  test  samples  when  subjected  to  monotonic  four- 
point  bending.  The  displacement  at  yield  was  1.0  mm  for 
the  samples  with  a  height  of  2.5  mm  and  1.07  mm  for  the 
samples  with  a  height  of  2.2  mm. 

A  mechanical  measure  of  the  amount  of  damage  that 
accumulated  within  the  test  sample  was  quantified  from  the 
magnitude  of  stiffness  degradation.  For  each  diagnostic 
cycle,  stiffness  was  calculated  from  a  linear  regression  of  the 
initial  portion  of  the  load-deformation  curve.  Specimen 
stiffness  decreased  nonuniformly  with  each  cycle  revealing 
increasing  amounts  of  damage  induced  within  each  cycle 
and  an  overall  damage  accumulation  by  the  end  of  the  pro¬ 
tocol  (Fig.  3B).  At  the  end  of  the  test  sequence,  the  overall 
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A. 


Damage  Parameter:  D  =  1  -  5^^  /  Sq  D  =  0  (no  damage) 

D  =  1  (fracture) 


B. 


DAMAGE  CYCLE 

FIG.  3.  (A)  A  15-cycle  loading  protocol  was  used  to  induce  dam¬ 

age  within  machined  cortical  bone  specimens  and  to  measure  re¬ 
sultant  stiffness  degradation.  Damage  was  induced  during  cycles 
2, 4,  6,  8, 10, 12,  and  14  by  conducting  relaxation  tests  at  increasing 
levels  of  applied  displacement  (expressed  as  a  percentage  of  the 
displacement  at  yield).  Diagnostic  cycles  were  interposed  between 
damaging  cycles.  Each  damage  cycle  was  preceded  by  a  diagnostic 
cycle  at  50%  of  displacement  at  yield.  A  5-minute  recovery  period 
was  introduced  after  the  damage  cycles  to  relieve  residual  internal 
stresses.  (B)  The  change  in  stiffness  calculated  between  sequential 
diagnostic  tests  was  plotted  vs.  cycle  number  for  the  damage  tests. 
The  dashed  curve  represents  a  specimen  showing  little  stiffness 
degradation  (i.e.,  little  damage  accumulation).  The  solid  line  rep¬ 
resents  a  specimen  showing  large  stiffness  degradation  (i.e.,  more 
damage  accumulation). 


damage  parameter,  D,  was  calculated  by  comparing  the 
stiffness  of  the  first  and  last  diagnostic  tests  such  that: 

D  =  \  -  5i5/5o,  (4) 

where  5^5  is  the  stiffness  of  the  last  diagnostic  cycle  and  Sq 
is  the  average  stiffness  of  the  first  two  diagnostic  cycles  (5i, 
^3)  and  the  first  damage  cycle  {S^}- 

Statistical  analysis 

All  data  were  regressed  against  age  using  linear  regres¬ 
sion  analysis  to  identify  the  properties  that  varied  signifi¬ 
cantly  with  age  (GraphPad  Prism,  San  Diego,  CA,  USA). 
To  determine  whether  bone  morphology  was  related  to  tis¬ 
sue  level  material  properties,  partial  correlation  coefficients 
were  determined  between  each  geometrical  parameter 
(e.g.,  Zap,  S)  and  each  tissue  level  mechanical  prop¬ 

erty  (modulus,  strength,  total  energy,  postyield  strain,  dam- 
ageability)  while  taking  age  into  consideration  (Minitab, 
State  College,  PA,  USA).^^^^ 

RESULTS 

The  sample  population  showed  broad  ranges  of  body 
size,  body  stature,  and  bone  morphology  values  (Table  1). 
Modulus  and  strength  showed  little  variation  among  indi¬ 
viduals  (CV  =  9.73%  and  4.62%,  respectively).  However, 


Table  1.  Variation  in  Properties  Among  Young  Adult 
Male  Tibias 


Property  Mean  ±  SD  Range 


Age  (years) 

Body  weight  (kg) 

Body  height  (cm) 

Body  mass  index  (kg/m^) 

Tibia  length  (cm) 

Cortical  area  (mm^) 

AP  width  (mm) 

ML  width  (mm) 

AP  moment  of  inertia  (mm**) 
ML  moment  of  inertia  (mm"*) 
Polar  moment  of  inertia  (mm"*) 
AP  section  modulus  (mm^) 

ML  section  modulus  (mm^) 

AP  slenderness  (1/mm^/kg) 

ML  slenderness  (f/mm^/kg) 
Modulus  (GPa) 

Strength  (MPa) 

Total  energy  (MPa) 

Post-yield  strain 
Damage  parameter  (D) 


32.9  ±  10.4 

17-46 

83.8  ±  25.1 

57.2-158.8 

177.7  ±  4.3 

170.2-182.9 

26.7  ±  8.0 

17.1-51.7 

38.1  ±1.9 

34.4-40.4 

355.9  ±  55.2 

268.3-511.1 

31.2  ±2.5 

28.4-36.9 

24.3  ±  2.3 

19.9-30.6 

34,390  ±  10,149 

19,466-60,698 

17,250  ±  5,945 

8,246-34,734 

51,640  ±15,886 

27,713-95,432 

3,279  ±  819 

1,925-5,172 

4,188  ±  907 

2,785-6,237 

9.9  ±  2.0 

6.5-13.4 

7.7  ±1.5 

5.4-10.2 

17.0  ±1.7 

13.4-19.0 

130.7  ±  6.1 

120.8-144.5 

3.3  ±  0.9 

2.3-5.6 

0.026  ±  0.006 

0.016-0.039 

0.165  ±  0.038 

0.103-0.253 

AP,  antero-posterior  direction;  ML,  medial-lateral  direction. 


postyield  strain  (CV  =  24.0%),  total  energy  (CV  = 
26.4%),  and  the  damage  parameter  (CV  =  23.0%)  all 
showed  large  variability  among  the  samples.  Morphological 
measures  such  as  AP  width,  section  modulus,  and  the  polar 
moment  of  inertia,  J  (Fig.  4),  increased  linearly  with  body 
weight  (R^  =  0.59,  p  <  0.003)  and  body  mass  index  (BMI; 

=  0.57,  p  <  0.004),  but  were  independent  of  body  height 
{R?  =  0.01,  p  <  0.7).  These  relationships  did  not  change 
when  the  body  weight  values  were  corrected  for  age  (data 
not  shown).  Body  height  was  uncorrelated  with  body 
weight  (R^  =  0.01,  p  <  0.8),  indicating  that  the  sample 
population  consisted  of  individuals  with  similar  heights  but 
widely  varying  body  weights. 

Significant  age-related  changes  were  observed  for  the  tis¬ 
sue-level  mechanical  properties  and  the  size  of  the  tibia.  A 
significant,  positive  correlation  was  observed  between  tibia 
slenderness  in  the  AP  {R?  =  0.31,  p  <  0.02)  and  ML  {R?  = 
0.24,  p  <  0.05)  directions  and  age.  However,  /^p,  /^l^ 

J  did  not  vary  with  age,  suggesting  that  the  variation  in 
slenderness  with  age  was  due  largely  to  higher  body  weight 
and  BMI  {R}  =  0.29-0.32,  p  <  0.03)  values  for  the  older 
individuals.  Although  tissue  modulus  did  not  vary  signifi¬ 
cantly  with  age,  tissue  strength  (i?^  =  0.53,  p  <  0.001),  post¬ 
yield  strain  {R^  =  0.44,  p  <  0.004),  and  total  energy  {R?  = 
0.32,  p  <  0.002)  were  significantly  lower  for  the  older  indi¬ 
viduals.  Furthermore,  a  significant,  negative  correlation  was 
observed  between  the  damage  parameter  and  age  {R^  = 
0.41,  p  <  0.006).  This  data  suggested  that,  whereas  the  tibia 
became  more  slender  relative  to  body  size  with  age,  the 
cortical  tissue  became  progressively  less  strong  and  less 
ductile  (i.e.,  more  brittle)  with  age. 

The  correlation  analysis  showed  significant  correlations 
between  tibial  morphology  and  the  mechanical  properties 
that  characterized  tissue  brittleness  and  damageability 
(Table  2).  The  relationships  among  tissue-level  mechanical 
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FIG.  4.  The  average  cross-sectional  polar  moment  of  inertia  of 
the  tibia  from  males  increased  with  (A)  body  weight  and  (B)  body 
mass  index  but  not  with  (C)  body  height.  Removing  the  outlier 
changed  the  values  to  0.24  (p  <  0.05)  for  A  and  0.20  {p  <  0.09) 
for  B. 


properties  and  cross-sectional  morphology  were  linear.  Post¬ 
yield  strain  and  total  energy  increased  significantly  with 
AP  width  (Figs.  5A  and  5B).  Modulus  decreased  with  /^p 
(p  <  0.07),  J  (p  <  0.08),  AP  section  modulus  (p  <  0.05),  and 
ML  section  modulus  (Figs.  5C-5F).  Tissue  damageability 
increased  with  tibia  slenderness  in  the  AP  (p  <  0.05;  Fig.  6) 
and  ML  (p  <  0.09)  directions.  These  correlations,  which 
were  independent  of  age,  indicated  that  a  narrower  bone 


was  comprised  of  tissue  that  failed  in  a  more  brittle  man¬ 
ner  and  accumulated  more  damage. 

DISCUSSION 

The  results  of  this  study  revealed  that  the  tissue-level 
mechanical  properties  of  cortical  bone  varied  with  the  size 
of  the  tibia.  Positive  correlations  were  observed  between 
measures  of  bone  size  (AP  width)  and  measures  of  tissue 
ductility  (postyield  strain,  total  energy),  and  negative  cor¬ 
relations  were  observed  between  bone  size  (moment  of  in¬ 
ertia,  section  modulus)  and  tissue  modulus.  Many  of  these 
correlations  were  significant.  The  lack  of  significant  corre¬ 
lation  with  all  measures  of  bone  size  can  be  attributed 
largely  to  the  complex  shape  of  the  tibia.  The  tibia  has  a 
triangular  cross-section  and,  consequently,  measures  of 
width  correlated  significantly  with  mechanically  relevant 
traits  like  cortical  area  and  moment  of  inertia  but  explained 
only  50-80%  of  the  variability  in  these  measures  (data  not 
shown).  These  correlations  would  be  greater  if  the  cross- 
section  had  a  circular  shape.  The  variability  in  these  corre¬ 
lations  was  sufficiently  large  that  neither  the  linear  (width) 
traits  nor  the  integrated  traits  like  area  and  moment  of 
inertia  correlated  significantly  with  a  particular  tissue-level 
mechanical  property  simultaneously.  Nevertheless,  the  data 
indicated  that  bones  with  smaller  width  were  comprised  of 
stiffer  and  less  ductile  (i.e.,  more  brittle)  material  compared 
with  larger,  more  robust  bones.  The  correlation  between 
tissue  ductility  and  bone  size  may  help  explain  why  male 
military  recruits^^’^’^^^  and  male  athletes^^^^  with  narrow 
bones  show  a  higher  ineidence  of  stress  fractures  compared 
with  individuals  with  wide  bones. 

The  development  of  the  slenderness  index^^^^  was  for  a 
“normal”  range  in  height  and  weight  and  is  probably  not 
useful  beyond  this  range.  However,  the  morphological 
variation  observed  in  our  sample  population  was  consistent 
with  that  reported  for  military  recruits^^’^^  and  runners,*^^^^ 
and  height  and  weight  were  consistent  with  recent  national 
averages. As  expected,  bone  size  varied  with  body 
weight,*^^^^  but  did  not  vary  with  height  (Fig.  Thus, 

narrow  bones  came  from  less  heavy  individuals  who  were  of 
similar  height  as  those  with  wide  tibias.  Weight  varied  more 
than  height  for  our  sample  population  similar  to  that  ob¬ 
served  for  the  aged-matched  national  data.  Furthermore, 
the  variability  in  weight,  specifically  inclusion  of  one  outlier 
(Fig.  4),  did  not  affect  the  results  (i.e.,  the  heaviest  person 
did  not  have  an  unusual  slenderness  value).  Thus,  the  bones 
used  in  this  study  seem  to  be  an  appropriate  size  relative  to 
body  type. 

The  variation  in  long  bone  slenderness  has  been  attrib¬ 
uted  to  genetic  and  environmental  factors  influencing 
growth  and  development^^^^  and  has  been  implicated  as  a 
risk  factor  for  osteoporotic  fracture.*^^®^  To  be  relevant  for 
military  recruits,  the  sample  population  should  have  ranged 
in  age  between  18  and  25  years.  However,  for  the  age  range 
in  this  study,  the  tissue-level  mechanical  properties  varied 
linearly  with  age  and  were  easily  corrected  using  a  linear 
regression  method.^^^^  Consequently,  the  correlation  analy¬ 
sis  presented  here  provides  relevant  insight  into  the  rela¬ 
tionship  observed  between  bone  size  and  stress  fracture  risk 
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Table  2.  Partial  Correlation  Coefficients  Taking  Age  Into  Consideration 


Modulus 

Strength 

Total  energy 

PY  strain 

Damageability 

Cortical  area 

-0.24 

0.03 

0.37 

0.40 

-0.23 

(0.35) 

(0.90) 

(0.14) 

(0.11) 

(0.37) 

AP  width 

-0.09 

-0.03 

0.57 

0.70 

-0.16 

(0.74) 

(0.90) 

(0.02) 

(0.002) 

(0.55) 

ML  width 

-0.32 

-0.22 

0.34 

0.41 

-0.19 

(0.21) 

(0.39) 

(0.18) 

(0.11) 

(0.45) 

AP  moment  of  inertia 

-0.45 

-0.10 

0.25 

0.32 

-0.27 

(0.07) 

(0.70) 

(0.34) 

(0.21) 

(0.29) 

ML  moment  of  inertia 

-0.39 

-0.18 

0.22 

0.28 

-0.25 

(0.12) 

(0.50) 

(0.39) 

(0.28) 

(0.32) 

Polar  moment  of  inertia 

-0.43 

-0.13 

0.24 

0.31 

-0.27 

(0.08) 

(0.61) 

(0.35) 

(0.23) 

(0.30) 

AP  section  modulus 

-0.50 

-0.13 

0.10 

0.14 

-0.30 

(0.04) 

(0.62) 

(0.70) 

(0.60) 

(0.25) 

ML  section  modulus 

-0.47 

-0.05 

0.18 

0.23 

-0.33 

(0.06) 

(0.84) 

(0.49) 

(0.38) 

(0.20) 

Ap  slenderness 

0.36 

0.10 

0.09 

0.03 

0.48 

(0.15) 

(0.69) 

(0.72) 

(0.92) 

(0.05) 

ML  slenderness 

0.23 

-0.01 

0.04 

-0.05 

0.41 

(0.39) 

(0.96) 

(0.89) 

(0.84) 

(0.10) 

Pearson  correlation  coefficients 

are  shown  with  p  values 

in  parentheses.  Significant  correlations  are 

shown  in  bold.  Abbreviations 

are  as  shown  in 

Table  1. 


for  young  adult  males.  Further  studies  are  needed  to  deter¬ 
mine  if  this  relationship  holds  over  a  wider  (older)  age 
range. 

The  data  provide  a  new  paradigm  that  may  explain  how 
variation  in  bone  slenderness  contributes  to  stress  fracture 
risk.  Individuals  with  narrow  tibias  were  previously  thought 
to  show  increased  fatigue  damage  during  intense  training 
because  the  smaller  bone  size  would  lead  to  an  overload 
situation  (i.e.,  higher  tissue  level  stresses). This  inter¬ 
pretation  was  based  on  the  assumption  that  tissue  mechani¬ 
cal  properties  did  not  vary  among  individuals.  However,  the 
current  results  indicated  that  tissue-level  mechanical  prop¬ 
erties  do  vary  among  individuals.  Specifically,  the  data  sug¬ 
gest  that  there  are  at  least  two  important  tissue-level  me¬ 
chanical  property  variations  that  need  to  be  considered  to 
understand  why  bone  size  is  a  risk  factor  for  stress  fractures. 
Narrower  tibias  were  comprised  of  tissue  that  was  more 
brittle  (low  total  energy)  and  was  prone  to  accumulate 
more  damage  compared  with  tissue  from  wider  tibia.  Hav¬ 
ing  tissue  that  is  more  or  less  damageable  may  be  inconse¬ 
quential  during  day-to-day  activities.  However,  tissue-level 
mechanical  properties  like  total  energy  and  ductility  be¬ 
come  particularly  important  in  defining  the  response  of 
bone  to  an  extreme  loading  condition,  such  as  that  expected 
during  military  training  or  during  a  fall.  Total  energy  de¬ 
fines  the  amount  of  energy  required  to  break  a  bone  (im¬ 
portant  during  a  fall)  and  ductility  and  damageability  define 
the  amount  of  damage  accumulated  under  overload  or  re¬ 
petitive  loading  (important  during  military  training).  Fur¬ 
thermore,  tissue  stresses  would  be  expected  to  remain 
higher  for  narrow  tibias  loaded  in  bending  or  torsion.  Mo¬ 
ment  of  inertia  is  related  to  the  external  diameter  raised  to 
the  fourth  power.  Because  whole  bone  stiffness  and 
strength  are  correlated  with  moment  of  inertia,^^^’^^^  a  bone 
with  a  large  external  diameter  should  also  show  large  over¬ 


all  stiffness  and  strength  values.  However,  the  -30%  varia¬ 
tion  in  tissue  modulus  (Table  1)  did  not  fully  compensate 
for  the  -100%  variation  in  the  moment  of  inertia  or  the 
section  modulus  (Table  1).^^^^  Thus,  in  situ  damage  accu¬ 
mulation  may  elicit  a  biological  response  (remodeling)  that, 
coupled  with  the  higher  tissue  stresses,  exacerbates  the  fa¬ 
tigue  process. Consequently,  individuals  with  narrow 
tibia  may  be  at  higher  risk  of  stress  fractures  because  of 
higher  in  vivo  tissue  stresses  (overloading)  coupled  with 
tissue  that  is  more  prone  to  accumulating  damage. 

The  data  may  also  help  explain  why  age  is  another  risk 
factor  for  stress  fractures. Bone  strength,  postyield 
strain,  and  total  energy  decreased  over  the  17-  to  46-year-old 
age  range.  This  was  consistent  with  previous  studies^^'’^^^^ 
and  indicated  that  cortical  bone  becomes  less  ductile  (i.e., 
more  brittle)  and  weaker  with  age  and  that  these  changes 
began  early  in  life.  This  age-related  decline  in  strength  and 
ductility  is  thought  to  be  a  result  of  increased  mineraliza¬ 
tion  and  remodeling.^^^’^"*^  Thus,  even  in  the  young  adult 
age  range,  the  amount  of  damage  accumulated  under  vig¬ 
orous  loading  regimens  would  be  expected  to  increase  with 
age.  This  variation  in  tissue  ductility  may  increase  the  sus¬ 
ceptibility  of  stress  fracture  risk  for  recruits  that  enter  into 
military  training  at  an  older  age. 

These  results,  and  those  of  others,  indicated  that  not 
all  cortical  tissue  was  constructed  in  the  same  manner. 
The  mechanical  properties  of  cortical  tissue  vary  with 
age,U‘’--32-35)  across  species,^^^^  among  bones  of  the  same 
individual,*-^®^  and  among  anatomical  sites  within  the  same 
bone.^^^’^®^  Here  we  showed  that  the  average  mechanical 
properties  of  cortical  tissue  also  varied  as  a  function  of  the 
overall  size  of  the  bone.  This  coupling  between  bone  mor¬ 
phology  and  tissue -level  mechanical  properties  has  been 
attributed  to  an  adaptive  response  of  bone.^^"*’^^’^®^  In  this 
study,  smaller  tibia  bone  size  was  coupled  with  an  increase 
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FIG.  5.  (A)  Postyield  strain  and  (B)  total 
energy  correlated  with  AP  width.  Modulus 
decreased  with  (C)  Zap,  (D)  J,  (E)  AP  sec¬ 
tion  modulus,  and  (F)  ML  section  modulus. 
Data  were  age-corrected  based  on  a  linear 
regression  method.^^^^ 
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FIG.  6.  Damageability  correlated  with  AP  slenderness  suggest¬ 
ing  that  tibias  that  were  more  slender  relative  to  body  size  and 
stature  were  comprised  of  tissue  that  accumulated  more  damage. 
Data  were  age-corrected  based  on  a  linear  regression  method.^^^ 

in  tissue  modulus.  The  goal  of  this  adaptive  response  is  to 
ensure  that  morphology  and  quality  together  meet  me¬ 
chanical  demands.  This  coupling  was  observed  when  com¬ 
paring  bones  subjected  to  widely  varying  mechanical  de¬ 
mands  from  different  species^^^^  and  has  also  been  used  to 
explain  the  maturation  of  bone  during  growth.^'*'^^^  Our 


current  results  suggested  that  this  coupling  might  also  exist 
for  a  particular  bone  (tibia)  within  the  same  species  (hu¬ 
man).  Additional  studies  are  needed  to  determine  if  similar 
relationships  between  morphology  and  quality  exist  for 
other  long  bones  (femur,  humerus,  radius). 

The  relationship  between  morphology  and  tissue-level 
mechanical  properties  observed  in  the  human  skeleton  was 
consistent  with  that  observed  for  the  mouse  skeleton.*-^"^^  In 
both  the  mouse  and  human  skeletons,  genetic  heterogene¬ 
ity  leads  to  variability  in  adult  bone  morphology  and  tissue 
level  mechanical  properties.  A  comparison  of  femurs  from 
A/J  and  C57/BL6  (B6)  inbred  strains  showed  that  A/J  fe¬ 
murs  were  more  slender  than  B6  as  a  result  of  the  two 
strains  having  similar  bone  lengths,  but  A/J  having  a  sig¬ 
nificantly  smaller  cross-sectional  size  and  shape.^^"*^  Despite 
the  difference  in  bone  size,  the  two  strains  showed  similar 
whole  bone  stiffness  values.  The  variability  in  bone  slen¬ 
derness  was  inversely  related  to  mineral  content,  suggesting 
that  bone  morphology  and  mineral  content  were  coordi- 
nately  regulated  so  whole  bone  stiffness  appropriately 
matched  the  mechanical  demands  imposed  by  weight  bear¬ 
ing.  However,  as  a  result  of  regulating  mineral  content  to 
match  bone  size,  A/J  femurs  failed  in  a  brittle  manner  and 
showed  poor  fatigue  properties.  In  the  human  skeleton, 
smaller  bones  were  stiffer  and  less  ductile.  Thus,  a  recipro¬ 
cal  relationship  was  observed  between  bone  stiffness  and 
ductility  for  both  skeleton  systems.  This  reciprocal  relation¬ 
ship  has  been  extensively  reported  for  cortical  bone,^^^^  and 
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it  is  thought  to  be  a  result  of  the  nature  of  the  compositional 
and  structural  factors  that  can  be  modulated  on  a  biological 
level.^"^^^^ Although  variation  in  mineral  content  may  have 
explained  the  differences  in  brittleness  for  the  mouse  skel¬ 
eton,  we  expect  that  the  human  skeleton  will  be  more  com¬ 
plex  and  that  the  variation  in  tissue-level  mechanical  prop¬ 
erties  will  be  a  consequence  of  variable  composition 
(mineral,  collagen,  water)  as  well  as  microarchitecture  (la¬ 
mellae,  osteon  size,  porosity). 

The  calculated  bending  modulus  and  strength  values, 
which  were  determined  from  machined  bone  samples  and 
were  thus  quantified  in  a  manner  that  was  independent  of 
bone  size,  were  consistent  with  bone  tensile  properties,^^°^ 
as  expected.  Test  samples  were  randomly  selected  to  obtain 
representative  mean  values  for  each  tibia  and  the  variation 
in  mechanical  properties  within  each  tibia  was  similar  to  the 
variability  observed  across  tibias.  Thus,  we  believe  that  the 
mean  values  reported  here  represent  the  generalized  tissue- 
level  mechanical  behavior  for  each  tibia. 

Compared  with  back-calculating  tissue-level  mechanical 
properties  from  whole  bone  failure  tests,  the  current 
method  of  measuring  tissue-level  mechanical  properties  di¬ 
rectly  from  machined  samples  provided  a  broader  range  of 
mechanical  properties  that  were  needed  to  better  under¬ 
stand  why  bone  size  is  a  risk  factor  for  stress  fractures.  The 
mechanical  properties  included  measures  of  ductility  (i.e., 
postyield  strain,  total  energy)  as  well  as  an  independent 
measure  of  damageability  (i.e.,  the  damage  parameter). 
These  properties  were  chosen  because  they  were  relevant 
for  understanding  the  material  response  of  bones  subjected 
to  the  vigorous,  repetitive  loading  associated  with  military 
training  and  running.  Postyield  strain  and  total  energy  rep¬ 
resent  measures  of  tissue  ductility  and  were  assessed  to 
discriminate  between  ductile  and  brittle  failure  modes.  Ma¬ 
terials  that  fail  in  a  brittle  manner  show  low  postyield  strain 
and  total  energy  values.  Variation  in  the  ductility  of  cortical 
bone  arises  from  differences  in  the  initiation,  accumulation, 
propagation,  and  coalescence  of  damage  in  the  form  of  mi¬ 
crocracks. Variation  in  the  damage  parameter  re¬ 
flected  differences  in  the  amount  of  damage  accumulated 
within  the  tissue  and/or  differences  in  the  way  damage  de¬ 
graded  tissue  stiffness.  The  damage  parameter  correlated 
negatively  with  postyield  strain  and  total  energy  (R^  = 
0.22-0.25,  p  <  0.05)  indicating  that  cortical  tissue  that  failed 
in  a  brittle  manner  also  tended  to  have  higher  tissue  dam¬ 
ageability  or,  accumulate  more  damage.  Although  the  ex 
vivo  bending  tests  do  not  necessarily  reflect  the  in  vivo 
loads  imposed  on  the  tibia,^^^’®°“^^^  the  bending  loads  were 
expected  to  induce  a  combination  of  tensile,  compressive, 
and  shear  damage^®^^  that  may  be  sufficiently  complex  to 
represent  a  generalized  variation  in  bone  quality  among 
human  tibias. 

The  results  of  this  study  provide  new  insight  into  why 
bone  size  is  a  risk  factor  for  stress  fractures.  Stress  fractures 
are  believed  to  be  a  consequence  of  excess  damage  accu¬ 
mulation  following  intense,  repetitive  activities.  Biological 
processes  that  attempt  to  repair  the  damage  may  further 
weaken  the  tissue  because  the  increased  resorption  results 
in  increased  tissue  porosity.^^"*^  However,  the  actual  contri¬ 
bution  of  biological  repair  processes  to  stress  fracture  risk 


remains  unclear.*^^^^  Damage,  in  the  form  of  microcracks,  is 
the  expected  sequelae  of  repetitive  loading  following  nor¬ 
mal,  daily  activities.^^®^  Intense  loading  conditions,  such  as 
those  associated  with  military  training  and  long  distance 
running,  are  expected  to  further  increase  in  situ  damage 
accumulation  and  degrade  tissue-level  mechanical  proper- 
ties.^^^^  Therefore,  under  extreme  loading  conditions  (e.g., 
military  training),  variation  in  bone  quality,  specifically  tis¬ 
sue  damageability,  may  be  a  contributing  factor  to  the  in¬ 
creased  risk  of  stress  fracture  for  individuals  with  more 
slender  bones.  The  current  data  suggested  that  bone  mor¬ 
phology  could  be  used  as  a  predictor  of  tissue  fragility  and 
stress  fracture  risk  in  the  absence  of  available  noninvasive 
imaging  techniques  that  accurately  measure  bone  damage- 
ability. 
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Genetic  Variation  in  Bone  Growth  Patterns  Defines  Adult  Mouse 

Bone  Fragility 


Christopher  Price, ^  Brad  C  Herman,^  Thomas  Lufkin,^  Haviva  M  Goldman,^  and  Karl  J  Jepsen^ 


ABSTRACT:  Femoral  morphology  and  composition  were  determined  for  three  inbred  mouse  strains  between 
ages  E18.5  and  1  year.  Genotype-specific  variation  in  postnatal,  pubertal,  and  postpubertal  growth  patterns 
and  mineral  accrual  explained  differences  in  adult  bone  trait  combinations  and  thus  bone  fragility. 

Introduction:  Fracture  risk  is  strongly  regulated  by  genetic  factors.  However,  this  regulation  is  generally 
considered  complex  and  polygenic.  Therefore,  the  development  of  effective  genetic-based  diagnostic  and 
treatment  tools  hinges  on  understanding  how  multiple  genes  and  multiple  cell  types  interact  to  create  me¬ 
chanically  functional  structures.  The  goal  of  this  study  was  to  connect  variability  in  whole  bone  mechanical 
function,  including  measures  of  fragility,  to  variability  in  the  biological  processes  underlying  skeletal  devel¬ 
opment.  We  accomplished  this  by  testing  for  variation  in  bone  morphology  and  composition  among  three 
inbred  mouse  strains  from  E18.5  to  1  year  of  age. 

Materials  and  Methods:  Mid-diaphyseal  cross-sectional  areas,  diameters,  moments  of  inertia,  and  ash  content 
were  determined  for  three  strains  of  mice  with  widely  differing  adult  whole  bone  femoral  mechanical  prop¬ 
erties  (A/J,  C57BL/6J,  and  C3H/HeJ)  at  E18.5  and  postnatal  days  1,  7, 14,  28,  56, 112, 182,  and  365  (n  =  5-15 
mice/strain/age). 

Results:  Significant  differences  in  the  magnitude  and  rate  of  change  in  morphological  and  compositional  bone 
traits  were  observed  among  the  three  strains  at  each  phase  of  growth,  including  prenatal,  postnatal,  pubertal, 
and  adult  ages.  These  genotype-specific  variations  in  growth  patterns  mathematically  determined  how  varia¬ 
tion  in  adult  bone  trait  combinations  and  mechanical  properties  arose.  Furthermore,  six  bone  traits  were 
identified  that  characterize  phenotypic  variability  in  femoral  growth.  These  include  (1)  bone  size  and  shape  at 
postnatal  day  1,  (2)  periosteal  and  (3)  endosteal  expansion  during  early  growth,  (4)  periosteal  expansion  and 
(5)  endosteal  contraction  in  later  growth,  and  (6)  ash  content.  These  results  show  that  genetic  variability  in 
adult  bone  traits  arises  from  variation  in  biological  processes  at  each  phase  of  growth. 

Conclusions:  Inbred  mice  achieve  different  combinations  of  adult  bone  traits  through  genotype-specific  regu¬ 
lation  of  bone  surface  activity,  growth  patterns,  and  whole  bone  mineral  accrual  throughout  femoral  devel¬ 
opment.  This  study  provides  a  systematic  approach,  which  can  be  applied  to  the  human  skeleton,  to  uncover 
genetic  control  mechanisms  influencing  bone  fragility. 
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INTRODUCTION 

Genetic  factors  are  generally  acknowledged  to  play 
an  important  role  in  the  determination  of  skeletal  fra¬ 
gility.  Clinically,  genetic  variability  in  bone  qual- 

ity,^^^  adult  bone  morphology,^^^  and  the  kinetics  of  bone 
loss*^"*’^^  all  contribute  to  the  variation  in  fracture  risk.*^^^ 
Therefore,  genetic-based  screening  tools  that  identify  an 
individual’s  skeletal  growth  potential  and/or  potential  for 
bone  loss  could  be  used  for  the  early  diagnosis  of  osteopo¬ 
rotic  fracture  risk  and  for  tailored  treatment  regimens 
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thereafter.^^^  However,  these  genetic  factors  are  not  com¬ 
pletely  understood.  Additionally,  the  processes  of  bone 
growth  and  bone  loss  are  genetically  complex,  such  that  the 
modeling  and  remodeling  of  complex  structures  involves 
the  coordinated  expression  of  many  genes  by  multiple  cell 
types  (osteoblasts,  osteoclasts,  and  osteocytes).^^^  There¬ 
fore,  a  primary  challenge  for  developing  effective  genetic- 
based  diagnostic  tools  for  fracture  risk  will  be  to  understand 
how  multiple  genes,  in  multiple  cell  types,  interact  through¬ 
out  life  to  create  and  maintain  structures  that  are  mechani¬ 
cally  and  biologically  functional.^^^ 

To  study  the  functional  connections  between  genotype 
and  skeletal  fragility,  we  adopted  a  top-down  approach  that 
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reveals  the  functional  connections  between  genotype  and 
skeletal  fragility/®^  Our  approach  uses  long  bones  from  in- 
bred  mouse  strains  to  relate  genetic  variability  in  whole 
bone  mechanical  properties  to  the  underlying  biological 
processes  regulating  bone  growth  and  loss  in  a  systematic 
and  hierarchical  manner.  Previous  research  using  inbred 
mice  has  shown  genotype-specific  variability  in  skeletal 
morphology/^*^^^^  mechanical  properties/^"^’^^^ 

bone  formation  rates, and  responsiveness  to  mechani¬ 
cal  stimulation.^^*'^®^  Previously,  we  showed  that  individual 
inbred  mouse  strains  exhibited  a  particular  combination  of 
cortical  area  (Ct.Ar),  polar  moment  of  inertia  (J^),  and 
mineral  content  (ash  content)  that  together  explained  66- 
88%  of  the  genetic  variability  in  adult  whole  bone  mechani¬ 
cal  properties.*^®^  These  results  indicated  that  the  combina¬ 
tion  of  transverse  bone  size,  shape,  and  material  quality  was 
responsible  for  the  particular  repertoire  of  stiffness, 
strength,  and  toughness  that  characterized  each  strain.  The 
next  step  in  our  systematic  approach  was  to  determine  how 
genetic  variability  in  these  three  physical  bone  traits  (Ct.Ar, 
Jo,  and  ash  content)  arises  during  the  growth  process. 

Bone  traits  are  influenced  by  both  genetic  and  environ¬ 
mental  factors  during  growth,  and  specific  combinations  of 
bone  traits  are  established  to  satisfy  the  mechanical  de¬ 
mands  associated  with  weight  bearing.*-^**^  Long  bone  di- 
aphyses  are  hollow  cylinders  whose  size  and  shape  can  be 
mathematically  determined  by  the  relative  amounts  of  bone 
apposition  and  resorption  on  the  periosteal  and  endosteal 
surfaces.^^^^^"*^  We  hypothesize  that  genetic  variability  in 
the  rates  of  periosteal  and  endosteal  expansion  throughout 
life  (i.e.,  growth  patterns)  defines  genotype-specific  combi¬ 
nations  of  adult  bone  traits.  To  test  this  hypothesis,  we 
compared  changes  in  femoral  morphology  (Ct.Ar  and  J^) 
and  composition  (ash  content)  across  development  (from 
embryonic  day  18.5  [E18.5]  to  1  year).  Transverse  growth 
was  studied  because  these  traits  contribute  predominantly 
to  the  resistance  of  in  vivo  loads. Three  genetically  dis¬ 
tinct  inbred  mouse  strains,  A/J  (A),  C57BL/6J  (B6),  and 
C3H/HeJ  (C3H),  with  widely  differing  adult  bone  trait 
magnitudes,  were  examined  to  better  understand  the  ex¬ 
pansive  growth  of  bone  surfaces  and  how  this  growth  leads 
to  unique  combinations  of  adult  morphological  (Ct.Ar, 
Tt.Ar,  Ma.Ar,  and  J^,)  and  compositional  (ash  content)  trait 
values. We  show  that  variability  in  skeletal  growth  pat¬ 
terns  reveals  the  biological  and  genetic  control  mechanisms 
that  define  adult  physical  bone  traits  and,  consequently, 
bone  fragility. 

MATERIALS  AND  METHODS 

Animals 

Female  (n  =  20/strain)  and  male  (n  =  10/strain)  A/J  (A), 
C57BL/6J  (B6),  and  C3H/HeJ  (C3H)  mice  were  purchased 
from  Jackson  Laboratory  (Bar  Harbor,  ME,  USA)  at  6-8 
weeks  of  age  and  used  to  establish  breeding  colonies.  Mice 
were  fed  a  standard  mouse  chow  (Purina  Laboratory  Chow 
5001;  Purina  Mills)  and  water  ad  libitum  and  kept  on  a  12-h 
lightidark  cycle.  Matings  were  monitored,  and  after  4  weeks 
of  weaning,  the  offspring  were  separated  by  sex  and  housed 


with  four  to  five  mice  per  cage.  Female  offspring  were 
killed  at  selected  time-points  for  analysis:  these  included 
E18.5,  postnatal  day  1, 7, 14, 28,  56, 112, 182,  and  365  (Table 
1).  For  the  E18.5  age  group,  pregnant  females  were  identi¬ 
fied  by  the  first  appearance  of  a  vaginal  plug  (designated  as 
day  E0.5).  The  pups  were  extracted  18  days  later  through 
caesarian  section  and  killed  by  decapitation.  All  mice  1  day 
of  age  and  older  were  killed  by  carbon  dioxide  asphyxiation 
and  decapitation.  The  Mount  Sinai  Institutional  Animal 
Care  and  Use  Committee  approved  all  procedures  for  the 
treatment  of  mice. 

Morphological  traits 

Right  femurs  were  harvested  and  embedded  in  poly¬ 
methylmethacrylate  as  described  previously.^^"^^  Transverse 
sections  of  the  femoral  diaphyses  (150  p,m  in  thickness) 
were  obtained  using  a  low-speed  diamond-coated  watering 
saw  (Buehler,  Lake  Bluff,  IL,  USA)  from  a  site  immedi¬ 
ately  distal  to  the  third  trochanter.  The  sections  were  af¬ 
fixed  to  acrylic  slides,  polished  to  l-|jim  finish,  and  stained 
using  von  Kossa’s  method.^^®^  The  sections  were  imaged 
with  a  digital  camera  (Sony  Exwave  HAD,  3CCD  Camera; 
Sony)  attached  to  a  visible  light  microscope  (Zeiss  Ax- 
ioplan2;  Zeiss).  The  pixel  resolution  of  this  optical  imaging 
system  was  2.1  p,m.  An  advanced  image  analysis  software 
package  (IMAQ  Vision  Builder  6.0;  National  Instruments 
Corp.,  Austin,  TX,  USA)  was  used  to  quantify  the  cross- 
sectional  areas  (total  area  [Tt.Ar],  cortical  area  [Ct.Ar],  and 
marrow  area  [Ma.Ar]),  cortical  thickness  (Ct.Th),  and  the 
polar  moment  of  inertia  (J^).  Young  bones  (ages  E18.5 
through  2  weeks  of  age)  had  a  predominantly  woven  ap¬ 
pearance  with  considerable  cortical  porosity.  Therefore,  at 
these  time -points,  morphometric  measurements  included 
all  porosities  within  the  cortical  shell  and  were  considered 
apparent  traits.  Three  cross-sections  were  analyzed  per 
bone,  and  the  values  were  averaged.  To  establish  rates  of 
surface  growth,  the  time  normalized  change  in  Tt.Ar  and 
Ma.Ar  during  the  early  phase  of  growth  (1-28  days  of  age) 
and  the  later  phase  of  growth  (28-112  days  of  age)  were 
obtained  from  the  following  equations  for  all  strains: 

ATt.Ari_28  =  [Tt.Ar28  -  avg(Tt.Ari)]/28  days 

(1) 

ATt.Ar28_ii2  =  [Tt.Arii2  -  avg(Tt.Ar28)]/84  days 

(2) 

AMa.Ar]^_28  =  [Ma.Arjs  -  avg(Ma.Ari)]/28  days 

(3) 

AMa. Arjyia  ArMax-112  =  [Ma.Ar]^i2  “  avg(Ma.Arf^j3  ^rMax)]/ 

(112  -  Ma.ArMax)  days  (4) 

The  subscript  refers  to  the  age(s)  over  which  the  rate  was 
determined.  Ma.ArMax  corresponds  to  the  age  at  which 
maximum  Ma.Ar  was  obtained  for  a  given  strain.  This  age 
was  used  as  the  initial  time-point  to  calculate  the  rate  of 
marrow  infilling  during  later  growth.  avg(Tt.Arj)  and  avg- 
(Ma.Atj)  refers  to  the  strain  average  Tt.Ar  and  Ma.Ar 
value  at  time-point  L 
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Table  1.  Body  Weight  and  Morphological  Traits  for  Femoral  Diaphyses  of  Three  Inbred  Mouse  Strains  Between 

Embryonic  Day  18.5  (E18.5)  and  1  Year  of  Age 


Cortical 

Total  area 

Marrow 

Average  cortical 

Strain  (age) 

N 

Weight  (g) 

area  (mm^) 

(mnd) 

area  (mm^) 

(mni*) 

thickness  (mm) 

A 

E18.5 

5 

0.11  ±  0.02* 

0.16  ±  0.02*’* 

0.05  ±  0.01*’* 

0.004  ±  0.001* 

0.10  ±  0.02 

Id 

15 

1.3  ±0.1 

0.13  ±  0.02* 

0.30  ±  0.05* 

0.17  ±  0.03 

0.010  ±  0.004* 

0.08  ±  0.01 

7d 

12 

2.8  ±  0.5* 

0.22  ±  0.02 

0.50  ±  0.05*’* 

0.28  ±  0.04*’* 

0.028  ±  0.005 

0.11  ±0.01* 

14d 

8 

6.0  ±  0.4 

0.23  ±  0.05* 

0.59  ±  0.08*’* 

0.35  ±  0.05*’* 

0.035  ±  0.011*’* 

0.10  ±  0.02* 

28d 

8 

10.8  ±1.6 

0.37  ±  0.03* 

0.82  ±  0.06*’* 

0.45  ±  0.04*’* 

0.076  ±  0.012*’* 

0.13  ±0.01* 

56d 

9 

16.4  ±  2.3 

0.51  ±  0.09* 

0.97  ±  0.10*’* 

0.46  ±  0.03* 

0.119  ±  0.027*’* 

0.18  ±  0.03 

112d 

13 

22.9  ±  3.0* 

0.74  ±  0.04* 

1.15  ±0.8*’* 

0.41  ±  0.05*’* 

0.189  ±  0.023*’* 

0.25  ±  0.02*’* 

182d 

12 

23.2  ±  1.6*-* 

0.79  ±  0.04* 

1.25  ±0.06*’* 

0.46  ±  0.03*’* 

0.224  ±  0.024*’* 

0.26  ±  0.02*’* 

365d 

B6 

E18.5 

9 

24.7  ±  2.4 

0.79  ±  0.06* 

1.28  ±0.11*’* 

0.49  ±  0.07* 

0.228  ±  0.037*’* 

0.25  ±  0.02* 

10 

_ 

0.17  ±  0.04*’* 

0.33  ±  0.04*’* 

0.16  ±  0.04*’* 

0.014  ±  0.004*’* 

0.10  ±  0.03 

Id 

15 

1.3  ±0.1 

0.16  ±  0.03* 

0.35  ±  0.03* 

0.19  ±  0.04 

0.014  ±  0.003* 

0.09  ±  0.02 

7d 

12 

3.1  ±0.3* 

0.19  ±  0.05 

0.61  ±  0.05* 

0.41  ±  0.05*’* 

0.032  ±  0.008 

0.08  ±  0.02* 

14d 

12 

6.4  ±  0.5* 

0.24  ±  0.04* 

0.90  ±  0.08*’* 

0.66  ±  0.06*’* 

0.061  ±  0.012* 

0.08  ±0.01*’* 

28d 

6 

10.9  ±1.0 

0.39  ±  0.05 

1.24  ±  0.09*’* 

0.85  ±  0.06*’* 

0.133  ±  0.023*’* 

0.11  ±0.01*’* 

56d 

11 

16.8  ±  0.8 

0.57  ±  0.03* 

1.47  ±0.07*’* 

0.90  ±  0.05*’* 

0.226  ±  0.022*’* 

0.16  ±0.01* 

112d 

10 

20.9  ±  0.6* 

0.75  ±  0.03* 

1.55  ±0.06*’* 

0.80  ±  0.04*’* 

0.297  ±  0.026* 

0.21  ±  0.03*’* 

182d 

9 

21.7  ±  0.9* 

0.76  ±  0.05* 

1.48  ±  0.09*’* 

0.72  ±  0.06*’* 

0.283  ±  0.029* 

0.20  ±  0.04*’* 

365d 

10 

25.2  ±1.8 

0.82  ±  0.06* 

1.82  ±0.10* 

1.00  ±  0.07*’* 

0.387  ±  0.044* 

0.20  ±  0.02* 

C3H 

E18.5 

8 

— 

0.13  ±  0.01* 

0.24  ±  0.03*’* 

0.11  ±  0.02*’* 

0.007  ±  0.001* 

0.10  ±0.01 

Id 

15 

1.4  ±0.1 

0.15  ±  0.03 

0.32  ±  0.03 

0.16  ±  0.02 

0.012  ±  0.002 

0.09  ±  0.02 

7d 

9 

4.4  ±  0.5*-* 

0.23  ±  0.06 

0.56  ±  0.07* 

0.33  ±  0.04*’* 

0.032  ±  0.010 

0.10  ±  0.02 

14d 

10 

5.7  ±  0.4* 

0.31  ±  0.03*’* 

0.77  ±0.11*’* 

0.46  ±  0.10*’* 

0.059  ±  0.015* 

0.12  ±0.01* 

28d 

8 

10.4  ±1.5 

0.43  ±  0.05* 

0.97  ±  0.7*’* 

0.54  ±  0.04*’* 

0.107  ±  0.017*’* 

0.14  ±0.01* 

56d 

10 

15.6  ±2.2 

0.67  ±  0.07*’* 

1.11  ±0.09*’* 

0.44  ±  0.04* 

0.174  ±  0.029*’* 

0.21  ±  0.04* 

112d 

11 

22.4  ±  1.3 

1.08  ±  0.04*’* 

1.41  ±0.06*’* 

0.34  ±  0.03*’* 

0.310  ±  0.025* 

0.34  ±  0.01*’* 

182d 

9 

20.9  ±  0.9* 

1.08  ±  0.05*’* 

1.35  ±  0.05*’* 

0.26  ±  0.04*’* 

0.285  ±  0.024* 

0.37  ±  0.03*’* 

365d 

7 

23.7  ±  1.8 

1.29  ±0.17*’* 

1.71  ±  0.23* 

0.42  ±  0.10* 

0.448  ±0.115* 

0.35  ±  0.04*’* 

*  Significantly  different  from  A,  ANOVA  p  <  0.05. 

^  Significantly  different  from  B6,  ANOVA  p  <  0.05. 

*  Significantly  different  from  C3H,  ANOVA  p  <  0.05. 


Compositional  traits 

Ash  content  was  assessed  for  the  left  femurs  of  all  mice  2 
weeks  of  age  and  older.  These  femurs  had  previously  been 
loaded  to  failure  in  four-point  bending  (data  not  included). 
The  diaphyseal  pieces  were  retrieved  and  cleaned  of  all  soft 
tissues  under  a  stereomicroscope.  The  hydrated,  dried,  and 
ashed  weights  were  determined  as  described  previously.*^®^ 
Water  content  was  defined  as  the  hydrated  weight  minus 
dried  weight  and  expressed  as  a  percentage  of  hydrated 
weight.  Ash  content  was  determined  as  the  ash  weight  nor¬ 
malized  for  hydrated  weight. 

Statistics 

All  results  are  expressed  as  mean  ±  SD.  Differences  in 
trait  values  among  the  three  genotypes  at  each  time -point 
were  determined  with  a  one-way  ANOVA  and  a  Tukey’s 
posthoc  test.  A  p  value  less  than  0.05  was  considered  sig¬ 
nificant  (GraphPad  Software,  San  Diego,  CA,  USA). 

RESULTS 

Overall  changes  in  femoral  diaphyseal  structure 

All  three  strains  underwent  qualitatively  similar  changes 
in  bone  architecture,  size,  and  shape  of  the  mid-diaphyseal 


cortex  as  a  function  of  age  (Fig.  1).  For  example,  each  strain 
exhibited  an  entirely  woven  or  porous  tissue  microstructure 
at  E18.5.  However,  by  14  days  of  age,  the  woven  tissue  was 
located  only  in  the  posterior-lateral  quadrant,  and  by  28 
days  of  age,  the  entire  mid-diaphysis  exhibited  a  solid,  com¬ 
pact,  lamellar  tissue  structure.  However,  despite  similar 
overall  growth  patterns,  significant  variation  in  bone  size, 
shape,  and  composition  were  observed  among  the  strains 
throughout  growth  (Table  1). 

Prenatal  bone  growth 

Significant  differences  in  total  area  (Tt.Ar),  cortical  area 
(Ct.Ar),  and  polar  moment  of  inertia  (J^)  were  observed 
as  early  as  E18.5  (Fig.  2A;  ANOVA,  p  <  0.05).  B6  mice 
showed  the  largest  Ct.Ar,  Tt.Ar,  and  J^,,  whereas  A  mice 
showed  the  smallest.  C3H  mice  more  closely  resembled 
A  mice  at  E18.5.  Between  E18.5  and  1  day  of  age  none  of 
the  morphological  traits  changed  in  B6  mice.  In  contrast, 
A  and  C3H  mice  showed  increases  in  all  three  traits,  so  that 
by  1  day  of  age  only  A  mice  showed  significantly  smaller 
Tt.Ar,  Ct.Ar,  and  values  compared  with  B6  (Fig.  2B; 
p  <  0.05). 
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FIG.  1.  Variation  in  femoral  cross-sectional  morphology  across  development  for  female  A,  B6,  and  C3H  mice.  Sections  represent 
equivalent  locations  within  the  femoral  diaphysis  taken  immediately  distal  to  the  third  trochanter.  Note  change  from  a  porous,  woven 
structure  to  a  compact,  lamellar  structure  between  PI  and  P28  days  of  age.  Scale  bar  shown  applies  to  all  images. 


A 


A 


AGE  (DAYS) 
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FIG.  2.  Variation  in  femoral  cortical  area  (Ct.Ar),  total  area 
(Tt.Ar),  and  polar  moment  of  inertia  (!„)  among  A,  B6,  and  C3H 
mice  at  (A)  embryonic  day  18.5  (E18.5)  and  (B)  postnatal  day  1. 
a,  b,  and  c  indicates  a  significant  difference  from  A,  B6,  and  C3H 
mice,  respectively  (ANOVA,  p  <  0.05). 


AGE  (DAYS) 

FIG.  3.  Variation  in  (A)  cortical  area  (Ct.Ar),  (B)  polar  moment 
of  inertia  (J,,),  and  (C)  ash  content  across  development  for  A,  B6, 
and  C3H  femurs.  Significant  differences  among  strains  in  Ct.Ar, 
Jq,  and  ash  content  are  indicated  in  Table  1. 


Postnatal  changes  in  Ct.Ar  and  Jg 

After  birth,  Ct.Ar  and  increased  nonlinearly  to  112 
days  of  age  (Figs.  3A  and  3B).  Significant  differences  in 
Ct.Ar  and  values  were  observed  among  the  strains 


throughout  development  (Table  1);  however,  these  differ¬ 
ences  became  more  prominent  at  56  days  of  age.  For  Ct.Ar, 
A  and  B6  mice  showed  similar  values  beyond  7  days  of  age, 
whereas  C3FI  mice  showed  larger  Ct.Ar  than  both  A  and 


GENETIC  VARIABILITY  IN  BONE  GROWTH 


1987 


A 


AGE  (DAYS) 

c 


1.25 


AGE  (DAYS) 


B 


D 


FIG.  4.  Variation  in  femoral  total  area 
(Tt.Ar)  and  marrow  area  (Ma.Ar)  among  A, 
B6,  and  C3H  female  mouse  strains.  (A) 
Changes  in  total  area  across  development. 
(B)  Rate  of  change  in  total  area  during  early 
growth  (P1-P28  days  of  age)  and  later 
growth  (28-112  days  of  age).  (C)  Changes  in 
marrow  area  across  development.  (D)  Rate 
of  change  in  marrow  area  during  early 
growth  (P1-P28  days  of  age)  and  growth  af¬ 
ter  maximum  marrow  area  was  achieved 
(Ma.ArMa^-112  days  of  age).  Significant  dif¬ 
ferences  among  strains  in  Tt.Ar  and  Ma.Ar 
are  indicated  in  Table  1.  a,  b,  and  c  indicates 
a  significant  difference  from  A,  B6,  and  C3H 
mice,  respectively,  at  the  indicated  time- 
point. 


B6  after  14  days  of  age  (ANOVA,  p  <  0.05).  For  J^,  B6  and 
C3H  mice  showed  similar  values  beyond  56  days  of  age, 
whereas  A  mice  showed  smaller  values  compared  with  B6 
and  C3H  beyond  14  days  of  age  (p  <  0.05).  These  differ¬ 
ences  in  physical  bone  traits  were  accompanied  by  nearly 
identical  age-related  increases  in  body  weight  among  all 
strains  (Table  1). 

Postnatal  expansion  of  the  periosteal  surface 

Tt.Ar  increased  nonuniformly  with  age  in  a  genotype- 
dependent  manner  (Fig.  4A).  Each  strain  achieved  80-90% 
of  adult  Tt.Ar  and  periosteal  diameter  by  28  days  of  age.  In 
contrast,  all  three  strains  achieved  only  40-50%  of  adult 
weight,  Ct.Ar,  and  by  this  age.  These  data  indicated  that 
the  periosteum  expanded  rapidly  for  all  strains  early  in  life. 
Flowever,  the  rate  of  expansion  at  the  periosteum  (i.e.,  in¬ 
crease  in  Tt.Ar  per  day)  during  early  growth,  and  the  size  of 
Tt.Ar  at  28  days  of  age  differed  significantly  among  the 
strains  (B6  >  C3H  >  A;  ANOVA,  p  <  0.05;  Fig.  4B;  Table 
1).  Between  28  and  112  days  of  age,  the  average  rate  of 
increase  in  Tt.Ar  was  reduced  by  -80%  compared  with  the 
rates  observed  between  PI  and  P28  days  of  age.  During  the 
later  phase  of  growth  (28-112  days),  the  rate  of  Tt.Ar  ex¬ 
pansion  was  greatest  in  C3H  mice,  whereas  A  and  B6  mice 
showed  similarly  lower  rates  of  expansion  {p  <  0.05).  Age- 
related  changes  in  periosteal  diameter  mirrored  the 
changes  in  Tt.Ar,  as  expected  (data  not  shown). 

Postnatal  expansion  of  the  endosteal  surface 

Like  Tt.Ar,  marrow  area  (Ma.Ar)  also  changed  nonuni¬ 
formly  with  age  and  in  a  genotype-dependent  manner  (Fig. 
4C).  All  mouse  strains  exhibited  similar  Ma.Ar  values  at 
birth  and  exhibited  rapid  postnatal  marrow  expansion  be¬ 
tween  PI  and  P28  days  of  age.  During  this  early  growth 
phase,  B6  femurs  exhibited  the  greatest  rate  of  endosteal 
expansion  (resorption),  followed  by  C3H  mice,  and  then  A 


mice,  resulting  in  significant  differences  in  Ma.Ar  size  at  28 
days  of  age  (A  <  B6  <  C3H;  ANOVA,  p  <  0.05;  Fig.  4D). 
Maximum  Ma.Ar  values  were  achieved  for  A  and  C3H 
mice  at  -28  days  of  age,  whereas  this  did  not  occur  until  56 
days  of  age  for  B6  mice.  After  reaching  maximal  Ma.Ar 
values,  all  strains  exhibited  decreases  in  Ma.Ar,  indicating 
that  the  endosteal  surface  underwent  a  reversal  in  biologi¬ 
cal  activity  from  net  resorption  to  net  apposition  (i.e.,  in¬ 
filling).  During  the  period  between  the  time  of  maximal 
Ma.Ar  and  112  days  of  age,  the  three  strains  of  mice 
showed  significantly  different  rates  of  endosteal  contraction 
(apposition;  A  <  B6  <  C3H;  p  <  0.05).  As  a  result,  adult 
Ma.Ar  at  112  days  of  age  was  largest  in  B6  mice  because  of 
the  early  rapid  expansion  of  the  marrow  cavity  and  smallest 
in  C3H  mice  because  of  its  slower  marrow  excavation  and 
larger  fractional  infilling.  Age-related  changes  in  endosteal 
diameter  were  similar  to  those  observed  for  Ma.Ar,  as  ex¬ 
pected  (data  not  shown). 

Age-related  changes  in  cortical  thickness 

Changes  in  cortical  thickness  (Ct.Th)  reflected  the  rela¬ 
tive  rate  of  expansion  of  the  periosteal  and  endosteal  sur¬ 
faces.  The  average  cortical  thickness  increased  nonlinearly 
between  birth  and  1  year  of  age  in  all  three  strains  (Fig.  5). 
Between  E18.5  and  7  days  of  age,  the  average  cortical  thick¬ 
ness  did  not  change  significantly  for  any  of  the  strains  in¬ 
dicating  that  expansion  of  the  marrow  cavity  kept  pace  with 
expansion  of  the  periosteum.  However,  beyond  14  days  of 
age,  cortical  thickness  increased  linearly  until  peak  values 
were  obtained  at  112  days  of  age.  After  14  days  of  age,  C3H 
femurs  exhibited  cortices  with  the  largest  average  thickness, 
A  femurs  exhibited  intermediate  cortical  thickness,  and  B6 
the  smallest  (ANOVA,  p  <  0.05). 

Postnatal  changes  in  ash  content 

In  all  three  strains,  ash  content  increased  with  age  to  112 
days  (Fig.  6).  Genotype-specific  variation  in  ash  content 
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FIG.  5.  Variation  in  femoral  cortical  thickness  (Ct.Th)  among  A, 
B6,  and  C3H  mice  across  development.  Significant  differences  in 
Ct.Th  among  strains  are  indicated  in  Table  1. 
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FIG.  6.  Variation  in  femoral  ash  content  across  development  for 
A,  B6,  and  C3H  femurs.  *B6  <  A  (ANOVA;  p  <  0.05);  **B6  < 
C3H  (ANOVA;  p  <  0.05). 


was  observed  at  the  earliest  time-point  assayed  (2  weeks  of 
age),  and  these  relationships  persisted  throughout  the  re¬ 
mainder  of  the  study. 


DISCUSSION 

This  study  showed  that  genetic  variability  in  inbred 
mouse  adult  bone  traits  (Ct.Ar,  J^,  and  ash  content)  results 
directly  from  genotype-specific  differences  in  periosteal  and 
endosteal  expansion  and  whole  bone  mineral  accrual,  as 
hypothesized.  Previous  studies  comparing  B6  and  C3H 
mouse  growth  suggested  that  genetic  variation  in  adult 
bone  mass  could  be  explained  by  prepubertal  and  pubertal 
growth  patterns^^^^  and  by  different  rates  of  bone  formation 
during  adolescent  growth.*^^®^  Additional  studies  examined 
B6  skeletal  development  aloneP^’^^^  However,  these  stud¬ 
ies  did  not  include  prenatal  and  immediate  postnatal 
growth  (<7  days).  Therefore,  it  remained  unclear  how  in 
utero  growth,  early  bone  cell  (osteoblast  and  osteoclast) 
activities,  and  peripubertal/adolescent  growth  contributed 
to  genetic  variability  in  adult  trait  magnitudes.  By  looking 
across  growth,  including  several  time-points  not  previously 
studied,  we  were  able  to  connect  early  and  later  bone  de¬ 
velopment.  Thus,  this  work  filled  in  gaps  from  previous 


studies  and  provided  new  insight  into  the  genetic  regulation 
of  adult  bone  traits.  Furthermore,  we  studied  both  the  spa¬ 
tial  distribution  of  bone  (J^)  and  matrix  mineralization  (ash 
content).  Previous  work  from  our  laboratory  has  shown 
that  these  traits  in  combination  with  tissue  amount  (Ct.Ar) 
are  important  in  the  determination  of  whole  bone  strength 
and  brittleness,  respectively.*^®^  Last,  we  included  a  third 
strain  (A),  which  was  previously  shown  to  have  a  small 
adult  bone  morphology  and  high  adult  ash  content,*®’*^"*^  to 
further  generalize  how  genetic  background  influences  adult 
bone  traits. 

Our  examination  of  bone  morphology  from  E18.5  to  365 
days  of  age  revealed  that  A,  B6,  and  C3H  inbred  mouse 
strains  showed  qualitatively  similar  transverse  femoral 
growth  patterns.  The  early  rapid  expansion  of  the  periosteal 
and  endosteal  surfaces  followed  by  marrow  infilling  ob¬ 
served  for  the  mouse  was  consistent  with  that  observed  for 
human  long  bones. As  such,  studying  the  genetic  regu¬ 
latory  pathways  controlling  murine  long  bone  growth  and 
the  attainment  of  peak  bone  properties  should  prove  valu¬ 
able  in  understanding  similar  processes  in  the  human  skel¬ 
eton. 

The  attainment  of  a  particular  combination  of  Ct.Ar  and 
Jo  in  adult  inbred  mice  can  be  explained  in  mathematical 
terms  based  on  variability  in  the  relative  expansion  of  the 
periosteal  and  endosteal  surfaces  during  growth  when  the 
femoral  diaphysis  is  modeled  as  a  simple  hollow  cylin- 
der/20-22)  gggajjsg  Ct.Ar  and  are  related  to  periosteal 
and  endosteal  diameters  differently, 

Ct.Ar  =  Tt.Ar  -  Ma.Ar  oc  Diap^^T  -  Dia^jj^jo^  (5) 

Jo  oc  Tt.Ar^  -  Ma.Ar^  cx  Diaporf"*  -  DiaEodo'^  (6) 

cylinders  of  different  diameters  will  necessarily  have  differ¬ 
ent  combinations  of  Ct.Ar  and  Jo  values.  Thus,  variability  in 
adult  Ct.Ar  and  Jo  arise  during  growth  necessarily  because 
differences  in  the  relative  expansion  of  the  periosteal  and 
endosteal  surface  (i.e.,  growth  patterns)  lead  to  different 
periosteal  and  endosteal  diameters. 

We  found  that  the  unique  pattern  of  periosteal  and  end¬ 
osteal  growth  exhibited  throughout  life  by  each  genotype 
could  indeed  explain,  in  mathematical  terms,  how  variabil¬ 
ity  in  adult  morphological  trait  combinations  (Ct.Ar  and  Jo) 
arises.  For  B6  mice,  the  large  Jo  and  small  Ct.Ar  observed 
in  adulthood  resulted  primarily  from  a  large  periosteal  and 
endosteal  expansion  during  early  growth  (before  puberty). 
For  C3H  mice,  an  intermediate  rate  of  periosteal  and  end¬ 
osteal  expansion  during  early  growth  (compared  with  A 
and  B6)  and  a  significantly  greater  periosteal  expansion  and 
marrow  infilling  after  puberty  led  to  the  large  Ct.Ar  and  Jo 
values  of  adult  C3H  mice.  As  a  result,  Ct.Ar  and  Jo  in¬ 
creased  at  a  faster  rate  in  C3H  mice  and  led  to  a  higher 
Ct.Ar  relative  to  Jo  because  a  larger  proportion  of  new 
bone  was  added,  in  a  less  mechanically  efficient  manner,  to 
the  endosteum.  It  is  noteworthy  that  B6  and  C3H  mice 
achieved  similar  adult  Jo  values.  However,  this  occurred 
through  very  different  growth  patterns  and,  presumably, 
different  biological  control  mechanisms.  For  A  mice,  the 
small  adult  Ct.Ar  and  Jo  were  a  result  of  having  the  smallest 
periosteal  and  endosteal  expansion  both  before  and  after 
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puberty.  Interestingly,  A  and  B6  mice  achieved  similar 
adult  Ct.Ar  values,  but  this  occurred  through  very  different 
growth  patterns,  suggesting  that  A  and  B6  mice  achieve 
similar  adult  Ct.Ar  through  different  biological  control 
mechanisms.  Together  these  results  show  that  inbred  mice 
build  morphologically  distinct  femurs  by  regulating  bone 
surface  activity  in  a  genotype-specific  manner. 

From  these  data,  we  identified  six  bone  traits  that  helped 
characterize  the  variability  in  transverse  femoral  develop¬ 
ment  among  three  strains  of  inbred  mice  and  that  provided 
additional  insight  into  the  biological  control  mechanisms 
underlying  variability  in  adult  bone  traits.  The  six  traits 
included  (1)  the  initial  size  of  the  bone  at  1  day  of  age,  (2) 
periosteal  expansion  during  the  early  (prepubertal),  rapid 
growth  phase,  (3)  endosteal  expansion  during  the  early 
growth  phase,  (4)  periosteal  expansion  during  the  later 
(peri-  to  postpubertal)  growth  phase,  (5)  marrow  infilling 
during  the  later  growth  phase,  and  (6)  ash  content  (material 
composition)  throughout  growth.  Interstrain  variability  in 
the  first  five  indices  reflects  differences  in  bone  cell  activity 
during  different  growth  phases,  and  the  last  index  reflects 
variation  in  the  tissue  mineralization  process.  Together, 
these  bone  traits  represent  phenotypic  indices  that  will  be 
useful  for  identifying  the  biological  (and  genetic)  mecha¬ 
nisms  that  influence  adult  whole  bone  mechanical  proper¬ 
ties. 

Variability  in  bone  size  and  shape  (Ct.Ar,  J^,  and  TAr)  at 
E18.5  and  1  day  of  age  show  that  these  inbred  mice  begin 
postnatal  skeletal  development  with  different  initial  mor¬ 
phologies.  These  results  suggest  that  bone  growth  and  adult 
bone  mass  might  be  related  to  allelic  variability  in  the  fac¬ 
tors  involved  in  in  utero  development.  Differences  in  tissue 
patterning  and  condensation  have  been  observed  in  the  cra- 
nio-facial  development  of  B6  and  C3H  mice,^^“^  and  similar 
variability  in  tissue  patterning  and  condensation,  cartilage 
formation,  or  initial  bone  formation  could  influence  long 
bone  morphology  at  birth.  Research  in  humans  has  identi¬ 
fied  correlations  between  in  utero  growth,  postnatal 
growth,  and  adult  bone  mass  that  support  the  idea  that 
pre-  and  postnatal  growth  influence  adult  skeletal  proper¬ 
ties. Additional  studies  in  mice  may  be  useful  in  es¬ 
tablishing  the  relationship  between  prenatal  morphology 
and  growth  and  adult  bone  trait  values. 

Investigation  of  early  growth  (1-28  days  of  age)  revealed 
that  prepubertal  growth  played  a  critical  role  in  defining 
adult  trait  values.  Strain  specific  differences  in  adult  exter¬ 
nal  size  were  established  during  early  growth,  suggesting 
that  variation  in  Tt.Ar,  bone  width,  and  maximal  Ma.Ar  are 
likely  the  result  of  allelic  variations  influencing  the  pro¬ 
cesses  controlling  rapid  surface  expansion  before  puberty. 
Measures  of  adult  bone  morphology,  which  have  been  cor¬ 
related  with  fracture  risk,^^’^"*^  are  also  largely  determined 
before  puberty  in  humans.^^^^  Therefore,  the  mouse  may 
provide  a  useful  model  for  studying  how  fracture  risk  is 
associated  with  the  biology  of  early  bone  growth  and  the 
establishment  of  external  size. 

Variability  in  ash  content  was  also  established  at  the  ear¬ 
liest  time-point  studied  (2  weeks  of  age)  and  was  consistent 
with  the  variability  observed  in  adulthood.  This  suggested 
that  differences  in  mineralization  and  material  composition 


among  the  strains  might  be  observed  even  earlier  in  life, 
possibly  before  birth,  or  within  the  first  few  days  of  post¬ 
natal  growth.  We  previously  showed  that  a  4.0%  variation 
in  ash  content  among  inbred  strains  was  correlated  with  a 
-3.4-  and  -1. 9-fold  variation  in  whole  bone  ductility  and 
toughness,  respectively. Therefore,  the  consistently 
lower  B6  femoral  ash  content,  compared  with  A  and  C3H 
mice,  implied  that  B6  femurs  were  constructed  of  bone  tis¬ 
sue  with  a  lower  material  stiffness  but  greater  ductility^^^^ 
throughout  life.  Together,  these  results  suggest  that  differ¬ 
ent  strains  of  inbred  mice  might  harbor  allelic  variations  in 
the  genetic  factors  inherent  to  or  governing  mineral  accrual. 

The  rate  and  extent  of  periosteal  and  endosteal  expan¬ 
sion  during  early  growth  appeared  to  be  well  matched  to 
ash  content  values,  suggesting  that  material  composition 
may  be  functionally  coupled  to  bone  growth  patterns.  B6 
mice,  with  the  most  rapid  prepubertal  surface  expansions 
and  largest  Tt.Ar  and  values  at  28  days  of  age,  showed 
small  ash  values.  In  contrast,  the  slower  expanding  bones  of 
A  and  C3H  mice  showed  larger  ash  values.  From  a  me¬ 
chanical  basis,  these  results  are  not  entirely  surprising  be¬ 
cause  the  addition  of  an  equivalent  amount  of  tissue  on  the 
outer  surface  of  a  large  bone  will  result  in  a  greater  increase 
in  Jq,  and  therefore  bone  stiffness,^^®'^^^  compared  with  a 
similar  addition  on  a  smaller  bone.  Thus,  the  rapid  expan¬ 
sion  of  the  periosteum  and  endosteum  in  B6  mice  estab¬ 
lished  a  more  efficient  growth  pattern  compared  with  A 
and  C3H  mice.  This  greater  structural  efficiency  would  sub¬ 
stantially  increase  whole  bone  stiffness  and  presumably  af¬ 
ford  B6  mice  the  ability  to  reach  a  mechanically  functional 
end  state  with  material  of  reduced  ash  content  and  thus 
reduced  tissue  level  stiffness. 

The  variability  in  femoral  growth  patterns  observed  after 
the  onset  of  puberty  were  entirely  consistent  with  previ¬ 
ously  reported  data.^^®’^^’^^^  Furthermore,  our  comparison 
of  early  and  later  growth  suggests  that  female  C3H  mice 
might  harbor  allelic  variations  that  significantly  alter  their 
growth  patterns  after  puberty  and  predispose  them  to  dras¬ 
tically  increased  bone  growth  in  association  with  sexual 
maturation.  Additional  studies  are  required  to  test  if  pu¬ 
bertal  maturation  and  sex  hormones  differentially  affect  the 
biological  processes  of  skeletal  growth  and  contribute  to  the 
variability  in  adult  bone  traits  among  these  strains. 

Significant  variation  in  the  aging  of  the  inbred  mouse 
skeleton  was  also  observed.  A  mice  showed  no  changes  in 
Tt.Ar  and  Ma.Ar  values  between  182  and  365  days  of  age, 
suggesting  little  periosteal  expansion  and  endosteal  expan¬ 
sion  during  this  time,  whereas  B6  and  C3H  mice  showed 
large  increases  in  Tt.Ar  and  Ma.Ar.  The  age-related  expan¬ 
sion  of  B6  and  C3H  femurs  is  similar  to  that  observed  for 
the  human  skeleton.*^^^’^®^  These  results  suggest  that  the 
biological  processes  associated  with  skeletal  aging  vary  with 
genetic  background. 

Although  female  A,  B6,  and  C3H  mice  exhibited  a  wide 
range  of  skeletal  phenotypes,  further  studies  are  warranted 
to  test  how  growth  patterns  vary  among  other  strains  of 
mice,  how  growth  patterns  vary  between  male  and  fe- 
males,*^^^’^^’^"*^  and  how  each  phase  of  growth  is  related  to 
biomechanical  properties.  In  a  standard  cylindrical  model 
of  growth,  measures  of  Tt.Ar  and  Ma.Ar  are  sufficient  to 
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describe  how  a  single  species  of  bone  cells  (osteoblasts  or 
osteoclast)  influence  surface  changes.  However,  the  trans¬ 
verse  growth  of  the  femoral  mid-diaphysis  occurs  through 
both  surface  expansion  and  cortical  drift.^^^’^^^  In  regions  of 
long  bones  exhibiting  drift,  both  apposition  and  resorption 
occur  at  different  locations  along  each  surface;  therefore, 
the  relative  contribution  of  osteoblasts  and  osteoclasts 
working  on  both  the  periosteum  and  endosteum  must  be 
considered  in  the  analysis  of  Tt. Ar  and  Ma.Ar,  respectively. 

From  this  study,  we  conclude  that  different  strains  of 
inbred  mice  build  distinct,  yet  mechanically  functional,  fe¬ 
murs  by  regulating  bone  surface  activity  and  mineral  ac¬ 
crual  in  a  genotype-specific  manner.  These  genotype- 
specific  growth  patterns  led  directly  to  particular 
combinations  of  adult  bone  traits,  and  thus,  specific  reper¬ 
toires  of  whole  bone  mechanical  properties.  By  looking 
across  all  ages  of  growth,  from  prenatal  development  to 
adulthood,  we  identified  six  bone  traits  that  can  be  used  to 
explain,  in  a  systematic  manner,  the  differences  in  femoral 
growth  patterns  among  inbred  mice.  This  characterization 
of  femoral  growth  for  A,  B6,  and  C3H  mice  throughout  life 
revealed  how  and  when  specific  combinations  of  adult  bone 
traits  arise  in  genetically  distinct  animals.  Thus,  the  inbred 
mouse  provides  a  valuable  model  to  determine  how  genetic 
background  influences  the  biological  processes  involved  in 
establishing  variability  in  adult  bone  traits,  and  thus,  whole 
bone  mechanical  properties. 
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Abstract 

Understanding  how  growth  influences  adult  bone  morphology  and  tissue-quality  should 
provide  important  insight  into  why  females  show  a  greater  incidence  of  stress  fractures  early  in 
life  and  fragility  fractures  later  in  life  compared  to  males.  The  objective  of  this  study  was  to  test 
whether  females  acquire  similar  tissue-level  mechanical  properties  as  males  by  the  time  peak 
bone  properties  are  established.  Standardized  beams  of  bone  were  machined  from  the  tibial 
diaphyses  of  14  young,  adult  females  ranging  in  age  from  22-46  years.  Data  for  males  (n=17, 
age=  17-46  years)  were  taken  from  a  prior  study.  Measures  of  tissue-level  mechanical  properties, 
including  stiffness,  strength,  ductility,  toughness,  and  damageability,  were  compared  between 
sexes  using  t-tests.  The  relationship  between  cross-sectional  morphology  and  tissue-level 
mechanical  properties  was  also  examined.  Males  and  females  showed  nearly  identical  tissue- 
level  mechanical  properties.  Both  sexes  also  showed  similar  age-related  degradation  of 
mechanical  properties  and  a  similar  relationship  between  cross-sectional  morphology  and  tissue- 
quality.  However,  for  all  body  sizes,  female  tibiae  were  smaller  relative  to  body  size  (i.e.,  less 
robust)  compared  to  males.  The  results  indicated  that  sex-specific  growth  patterns  affected 
transverse  bone  size,  but  did  not  affect  tissue  level  mechanical  properties.  This,  combined  with 
the  observation  that  young,  adult  female  long  bones  are  undersized  relative  to  body  size,  suggests 
that  adult  females  would  be  expected  to  accumulate  more  damage  under  intense  loading 
compared  to  males.  This  may  be  a  contributing  factor  to  the  greater  incidence  of  stress  fractures 
observed  for  female  military  recruits. 
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Introduction 

Understanding  why  females  show  a  greater  ineidenee  of  stress  fraetures  early  in  life  [1,2] 
and  fragility  fraetures  later  in  life  [3-5]  eompared  to  males  may  provide  important  insight  into 
the  biologieal  [6]  and  meehanieal  [7]  mechanisms  underlying  fracture  risk.  Sexual  dimorphism  is 
one  factor  that  may  contribute  to  this  discrepancy  [6,  8],  because  bone  size  and  shape  are  known 
determinants  of  bone  strength  [9].  Importantly,  fracture  also  depends  on  tissue-level  mechanical 
properties.  Although  it  is  known  that  sex-specific  growth  patterns  influence  bone  structure,  it  is 
unclear  whether  variable  growth  patterns  affect  the  construction  of  bone  matrix  in  a  way  that 
leads  to  differences  in  tissue-level  mechanical  properties.  The  few  studies  that  examined  tissue- 
level  mechanical  properties  of  long  bones  from  adult  males  and  females  never  directly  tested  for 
sex-specific  differences  [10,  11].  The  objective  of  this  study  was  to  test  whether  females  acquire 
similar  tissue-level  mechanical  properties  as  males  by  the  time  peak  bone  properties  are 
established.  In  addition  to  tissue-level  stiffness  and  strength,  we  also  quantified  tissue- 
damageability,  toughness,  and  ductility,  because  these  latter  mechanical  properties  provide 
insight  into  the  material  response  of  bone  when  subjected  to  extreme  load  conditions  such  as 
during  a  fall,  as  well  as  the  amount  of  damage  accumulated  within  the  bone  when  subjected  to 
intense  physical  activity.  We  also  tested  whether  females  show  a  correlation  between  cross- 
sectional  morphology  and  tissue-level  mechanical  properties  similar  to  that  previously  observed 
for  males  [12]. 

Materials  and  methods 

Sample  population 

Tibiae  of  14  female  donors  (12  Caucasian,  1  African-American,  1  unknown)  aged  22  -  46 
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years  (average  age  =  36.9  +  8.1  years)  were  acquired  from  the  Musculoskeletal  Transplant 
Foundation  (Edison,  NJ,  USA)  and  the  National  Disease  Research  Interchange  (Philadelphia, 
PA,  USA).  Donor  body  weight  and  height  were  obtained  from  the  source.  Only  tibiae  from 
donors  with  no  known  skeletal  pathology  were  included  in  the  study.  The  tibiae  were  freshly 
harvested,  wrapped  in  wet  gauze,  and  stored  in  plastic  bags  at  -40°C.  Data  for  male  tibiae  (n=17, 
average  age  =  32.9  +  10.4,  age  range  =  17-46  years)  were  taken  from  a  previous  publication  [12]. 

Whole  bone  morphology 

Tibia  length  (L)  was  measured  as  the  average  distance  between  the  middle  of  the  talar 
trochlear  facet  and  the  medial  and  lateral  proximal  condyles[13]  using  a  large-capacity  slide 
caliper  with  an  accuracy  of  +2. 5mm  (Mantex  Precision,  Haglof  Inc.,  Madison,  MS,  USA).  Tibia 
widths  in  the  antero-posterior  (WidthAp)  and  medial-lateral  (WidthMr)  directions  were  measured 
using  a  300mm  vernier  caliper  with  an  accuracy  of  +0.02mm  (Fowler  Company  Inc.,  Newton, 
MA,  USA).  Measurements  were  taken  at  10%  intervals  from  30%  to  70%  of  the  total  tibia  length 
and  averaged  in  order  to  assess  the  diaphyseal  morphological  traits  [14]. 

Cross-sectional  morphological  traits  were  quantified  from  3-mm  thick  mid-diaphyseal  cross- 
sections  cut  at  30,  50,  and  70%  of  the  total  tibia  length  using  a  diamond  coated  metallurgical  saw 
(Model  660,  South  Bay  Technology  Inc.,  San  Clemente,  CA,  USA).  Each  cross-section  was 
digitally  imaged  (0.024mm/pixel)  and  analyzed  using  IMAQ  Vision  Builder  (version  6.0, 
National  Instruments  Corp.,  Austin,  TX,  USA).  Morphological  traits  included  cortical  area 
(CtAr),  the  moments  of  inertia  about  the  antero-posterior  (Iap)  and  medial-lateral  (1ml)  axes,  the 
polar  moment  of  inertia  (J=  Iap  +  Iml),  and  the  section  modulus  in  the  AP  (J/WidthAp/2)  and  ML 
(J/WidthMiy2)  directions.  Moment  of  inertia  and  section  modulus  were  assessed  because  these 
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geometric  measures  are  related  to  the  bending  and  torsional  stiffness  of  intact  tibiae.  A 
slenderness  index  (S)  was  calculated  as  the  ratio  of  the  AP  and  ML  section  modulus  values, 
respectively,  to  tibia  length  and  body  weight; 

Sap  =  (l/J/(WidthAp/2))/(L*BW)  (1) 

Sml  =  (l/J/(WidthML/2))/(L*BW)  (2) 

where  L  =  tibia  length  (mm)  and  BW  =  body  weight  (kg).  For  semantic  reasons,  the  formulae  are 
the  inverse  to  that  used  previously  [15],  such  that  an  increase  in  Sap  or  Sml  indicates  a  more 
slender  (gracile)  bone.  All  morphological  traits  were  averaged  over  the  three  cross-sections  for 
each  tibia. 

Bone  sample  generation 

Cortical  bone  samples  were  prepared  for  biomechanical  testing  by  machining  beams  with  a 
rectangular  cross-section  from  the  diaphysis  of  each  tibia  using  an  automated  CNC  milling 
machine  under  constant  irrigation  (Modela  MDX-20,  Roland  DGA  Corp.,  Irvine,  CA,  USA). 
Sample  width  (circumferential  direction  of  tibia)  was  machined  to  5mm  and  length  (longitudinal 
direction)  was  machined  to  55mm  for  all  samples.  Sample  height  (radial  direction)  was  either  2.2 
or  2.5mm,  depending  on  the  cortical  thickness  of  the  tibia.  A  total  of  eight  samples  were 
generated  from  each  tibia  and  randomly  distributed  to  monotonic  (n  =  4)  and  damage 
accumulation  (n  =  4)  test  groups.  All  samples  were  stored  at  -40°C  in  gauze  saturated  with 
phosphate  buffered  saline  (PBS)  with  added  calcium  [16]  and  placed  individually  in  airtight 
bags. 
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Monotonic  failure  properties 

Tissue-level  meehanieal  properties  were  assessed  by  loading  four  eortieal  bone  samples  from 
eaeh  tibia  to  failure  in  4-point  bending  at  0.05mm/s  using  a  servohydraulie  materials  testing 
system  (Instron  model  8872,  Instron  Corp.,  Canton,  MA,  USA).  Speeimens  were  submerged  in  a 
PBS  solution  with  added  ealeium  and  maintained  at  37°C  throughout  all  tests.  Load  and 
defleetion  were  eonverted  to  stress  and  strain  as  deseribed  previously  [12].  These  bending 
equations  take  yielding  into  eonsideration  [17],  and  thus  provide  an  estimate  of  tissue  strength 
that  is  eonsistent  with  tensile  meehanieal  properties  [18].  All  properties  were  averaged  over  the 
four  samples  tested  for  eaeh  tibia. 

Damageability 

Tissue  damageability  was  assessed  using  the  same  protoeol  deseribed  previously  [12].  This 
protoeol  was  designed  to  induee  and  aeeumulate  eraeks  in  eortieal  bone  speeimens  and  to 
measure  the  degradation  of  stiffness  and  relaxation  as  surrogate  measures  of  damage 
aeeumulation  [19].  Four  eortieal  bone  samples  from  eaeh  tibia  were  subjeeted  to  a  fifteen-eyele 
damage  aeeumulation  protoeol.  For  this  protoeol,  diagnostie  eyeles  were  interposed  between 
damage  eyeles  to  assess  damage-assoeiated  degradation  at  eaeh  load  level.  For  the  diagnostie 
eyeles,  the  speeimens  were  loaded  in  four-point  bending  at  0.5mm/s  to  50%  of  the  average 
displaeement  at  yield  (determined  from  the  monotonie  tests),  held  for  60  seeonds,  and  then 
unloaded  at  0.5mm/s.  For  the  damage  eyeles,  the  speeimens  were  loaded  at  0.5mm/s  to  50,  75, 
100,  125,  150,  175,  and  200%  of  displaeement  at  yield  respeetively,  held  for  60  seeonds,  and 
then  unloaded  at  0.5mm/see.  A  5-minute  reeovery  period  followed  eaeh  damage  eyele. 
Displaeement  at  yield  was  used  as  a  referenee  in  the  damage  eyeles  beeause  this  parameter 
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showed  little  variation  among  the  test  samples  when  subjected  to  mono  tonic  four-point  bending. 
The  displacement  at  yield  was  1.0mm  for  the  samples  with  a  height  of  2.5mm  and  1.07mm  for 
the  samples  with  a  height  of  2.2  mm. 

Damage  accumulation  within  the  test  sample  was  quantified  by  calculating  the  stiffness 
degradation  and  the  change  in  relaxation.  We  included  measures  of  relaxation  in  this  analysis 
because  prior  research  showed  that  damage  accumulation  affects  viscoelastic  behavior  more 
dramatically  compared  to  stiffness  [19-22]  and  thus  would  provide  a  more  sensitive  assay  to 
assess  variation  in  tissue  damageability.  For  each  diagnostic  cycle,  i,  stiffness  {S^  was  calculated 
from  a  linear  regression  of  the  initial  portion  of  the  load-deformation  curve  and  relaxation  {R,) 
was  defined  as  the  decrease  in  load  over  the  60-second  hold  period.  The  total  amount  of  damage 
induced  in  each  sample  was  calculated  by  comparing  the  stiffness  and  relaxation  values 
measured  at  the  end  of  the  entire  test  sequence  {Sis,  Ris)  to  those  before  {So,  Ro).  To  test  if  the 
damage  accumulation  process  was  altered,  degradation  measures  were  also  calculated  after  each 
damage  cycle.  The  general  equations  for  degradation  measures  are  as  follows: 


Dstiff,-  =  1  —  SJ  So 

(4) 

DrELAX/  =  RJRo 

(5) 

where  Si  and  Rt  are  the  stiffness  and  total  relaxation,  respectively,  of  the  cycle  i,  and  So  and  Ro 
are  the  average  stiffness  and  total  relaxation,  respectively,  of  the  first  two  diagnostic  cycles  and 
the  first  damage  cycle. 

Statistical  analysis 

Mechanical  properties  and  morphological  traits  that  changed  significantly  with  age  were 
corrected  by  regression  analysis  [23].  Differences  between  females  and  males  were  determined 
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using  a  Student’s  t-test  with  corrections  applied  for  unequal  variances  (GraphPad  Software,  Inc; 
San  Diego,  CA  USA).  To  determine  if  female  and  male  tibiae  were  similarly  adapted  to  body 
size,  the  section  modulus  (cross-sectional  moment  of  inertia/bone  width)  was  regressed  against 
the  product  of  body  weight  and  tibia  length  [13]  and  the  slopes  were  compared  between  sexes 
(ANCOVA).  To  determine  if  males  and  females  show  a  similar  relationship  between  tibial  cross- 
sectional  morphology  and  tissue-level  mechanical  properties  [12],  linear  regressions  were 
performed  and  differences  in  slopes  and  intercepts  between  female  and  male  regressions  were 
determined  using  a  method  equivalent  to  an  Analysis  of  Covariance  (GraphPad  Software,  Inc; 
San  Diego,  CA  USA).  Since  no  differences  were  found  between  male  and  female  regressions,  a 
matrix  of  Pearson  correlation  coefficients  was  constructed  using  the  combined  female  and  male 
datasets  to  identify  the  morphology-tissue  quality  relationships  that  were  significant. 

Results 

Sample  Population 

Female  donors  were  shorter  than  males  (163  +  6  cm  vs.  178  +  4cm;  p<0.0001),  but  were  not 
different  in  body  weight  (74  +  21kg  vs.  84  +  25kg;  p<0.2). 

Monotonic  properties 

Females  and  males  showed  similar  regressions  for  tissue  modulus,  strength,  post-yield  strain, 
and  energy-to-failure  versus  age  (Fig.  1).  No  differences  in  the  average  tissue-level  mechanical 
properties  were  observed  between  females  and  males  (Table  1),  even  when  corrected  for 
differences  in  donor  age.  This  included  measures  of  stiffness  (modulus),  strength,  ductility 
(failure  strain,  post-yield  strain),  and  toughness  (energy-to-failure).  Thus,  males  and  females 
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showed  similar  tissue-level  monotonic  properties. 

Damageability 

No  differences  in  either  total  stiffness  degradation  or  the  total  change  in  relaxation  (Table  2) 
were  observed  between  females  and  males.  Further,  the  changes  in  stiffness  and  relaxation  after 
each  damage  cycle  were  similar  for  males  and  females  (Fig.  2).  These  results  indicated  that  the 
loading  protocol  induced  similar  amounts  of  damage  at  each  load  step  for  females  and  males. 

Morphology-Tissue  Quality  Relationships 

None  of  the  morphological  traits  showed  a  significant  correlation  with  age  for  either  males  or 
females  (data  not  shown).  Adult  females  achieved  a  different  bone  size  and  shape  compared  to 
males:  female  tibiae  were  smaller  (J,  CtAr,  Width),  more  slender  (Sap,  Sml,  J/L),  and  showed  a 
smaller  J:A  ratio  compared  to  males  (Table  3).  Females  and  males  showed  similar  slopes  for  the 
regression  of  J/WidthAp  versus  the  product  of  body  weight  and  tibial  length  (p<0.17,  ANCOVA). 
However,  females  showed  a  significantly  smaller  intercept  compared  to  males  (p<0.0001, 
ANCOVA),  indicating  that,  for  all  body  sizes,  female  tibiae  were  smaller  relative  to  body  size 
(i.e.,  less  robust)  compared  to  males  (Fig.  3).  Similar  results  were  observed  when  using 
J/WidthML. 

The  tissue-level  mechanical  properties  were  regressed  against  bone  morphological  traits  and 
no  significant  differences  in  slopes  were  found  for  females  and  males.  A  correlation  analysis 
(Table  4),  which  was  conducted  using  the  combined  male  and  female  data  sets,  showed  that 
tissue-level  modulus  decreased  with  increasing  bone  section  modulus  (J/WidthAp,  J/WidthMr), 
post-yield  strain  increased  with  increasing  bone  width  (WidthAp),  and  tissue-damage  (Dstiff) 
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increased  with  inereasing  tibial  slenderness  (Sap,  Sml)-  The  signifieant  regressions  are  shown  in 
Figure  4.  The  amount  of  overlap  between  the  female  and  male  datasets  depended  on  whether  the 
bone  morphology  traits  were  adjusted  for  body  size.  For  unadjusted  traits  like  width  and  J/width, 
the  female  data  simply  extended  the  relationships  observed  for  the  males  (Fig.  4A,B)-  For  the 
regression  involving  bone  slenderness,  whieh  is  adjusted  for  body  weight  and  tibial  length,  the 
female  data  overlapped  substantially  with  the  male  data  (Fig.  4C). 

Discussion 

Tissue-level  meehanical  properties,  ineluding  measures  of  fragility  and  damageability,  were 
assessed  for  eortieal  bone  samples  that  were  maehined  from  the  diaphyses  of  young,  adult  male 
and  female  tibiae.  The  data  indieated  that  the  tibial  diaphyses  of  females  and  males  were 
eomposed  of  material  having  similar  tissue-level  mechanical  properties,  and  that  the  meehanical 
properties  degraded  with  age  at  similar  rates  over  the  age  range  of  17-46  years.  Further,  males 
and  females  showed  a  similar  relationship  between  tissue-level  meehanieal  properties  and  eross- 
seetional  morphology.  These  data  suggested  that  the  genetic  [24-26]  and  environmental  [27] 
faetors  eontributing  to  sex-speeifie  growth  patterns  affeeted  adult  tibial  eross-seetional  size  and 
shape,  but  did  not  affect  bone  matrix  eonstruetion,  mineralization,  and  organization  in  a  way  that 
signifieantly  affeeted  tissue-level  meehanieal  properties.  This  analysis  ineluded  meehanieal 
properties  like  stiffness,  which  is  important  for  day-to-day  activities,  as  well  strength,  ductility, 
toughness,  and  damageability,  whieh  provide  insight  into  the  amount  of  damage  aeeumulated 
during  intense  physieal  aetivity  and  the  fraeture  resistanee  during  an  extreme  load  eondition  sueh 
as  a  fall. 

Our  data  eonfirm  that  males  and  females  achieve  similar  adult  tissue-level  meehanical 
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properties,  which  has  only  been  assumed  in  prior  analyses  [13].  Prior  studies  reported  property 
values  for  this  age-range,  but  did  not  directly  compare  data  for  females  and  males  [10,  11].  The 
mechanical  properties  measured  from  the  4-point  bending  tests  in  this  study  were  consistent  with 
the  mechanical  properties  of  cortical  bone  measured  previously  in  tension  [28-30],  which  was 
expected  given  that  failure  in  our  bending  tests  occurred  following  fracture  initiated  within  the 
tensile  region.  Further,  the  bending  formula  that  was  used  to  convert  load  and  deflection  to  stress 
and  strain,  because  it  accounted  for  nonlinear  effects  [17],  provided  a  more  appropriate  estimate 
of  bending  strength  compared  to  standard  beam  theory  [18].  Because  the  bone  samples  used  in 
this  study  were  machined  from  within  the  middle  of  the  cortex,  the  tissue  that  was  added  to  the 
subperiosteal  or  subendosteal  surfaces  during  and  after  puberty  was  likely  not  tested  in  our  study. 
Thus,  the  current  data  could  not  be  used  to  determine  if  sex-specific  differences  in  tissue-quality 
exist  following  puberty.  This  data  thus  reflects  the  tissue  properties  established  early  in  life,  plus 
the  modifications  to  these  properties  associated  with  osteonal-remodeling  and  cortical  drift  [31] 
during  the  ensuing  years. 

Although  the  current  data  were  collected  by  cross-sectional  means,  the  mechanical  property 
versus  age  regressions  implied  that  in  situ  tissue-level  mechanical  properties  begin  to  decline 
early  in  adulthood  for  males  and  females,  even  before  bone  loss  becomes  measurable  [32]. 
Similar  regressions  were  observed  previously  for  data  that  included  a  much  broader  age  range 
(19-90  years)  [10,  29,  30,  33].  Stiffness  and  strength  typically  increase  until  about  30  years  of 
age  and  then  decline  slowly  thereafter  [10].  In  contrast,  the  post-yield  behavior  of  bone  degrades 
rapidly  with  age,  and  this  loss  in  ductility  is  a  primary  reason  why  bone  tissue  becomes 
progressively  more  brittle  [30]  and  shows  lower  impact  strength  [11].  The  decline  in  tissue-level 
strength,  ductility,  and  damageability  have  been  attributed  to  age-related  changes  in  ash  content 
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[11] ,  osteonal  remodeling  [34],  collagen  quality  [35,  36],  mineral  density  distribution  [37],  and 
porosity  [11,  30].  Beyond  50  years  of  age,  the  mechanical  properties  of  females  and  males  may 
follow  different  degradation  pathways  associated  with  sex-specific  bone  loss  patterns  [30]. 

The  similarity  in  tissue-level  mechanical  properties  for  adult  men  and  women  may  be 
explained  based  on  early  transverse  bone-growth  patterns.  During  pre -pubertal  growth,  the 
periosteal  and  endosteal  surfaces  of  female  and  male  long  bone  diaphyses  follow  nearly  identical 
expansion  rates  [38].  Despite  the  dimorphic  transverse  bone-growth  patterns  arising  during 
puberty  [38],  the  long  bones  of  both  sexes  appear  to  be  constructed  in  a  way  so  that  by  20  years 
of  age  bone  size  is  properly  adapted  to  satisfy  loading  demands  associated  with  body  weight  [39- 
42].  Whole  bone  stiffness  and  peak  tissue  strains  are  thus  kept  at  proper  levels  [43]  during 
growth  simply  by  depositing  enough  tissue  on  the  periosteal  and  endosteal  surfaces  to  keep  pace 
with  weight  gain.  Thus,  because  there  is  no  discrepancy  in  the  relationship  between  bone  size 
and  body  size  during  growth,  males  and  females  can  construct  the  tibia  with  similar  tissue-level 
mechanical  properties  without  a  loss  in  mechanical  function. 

Although  each  of  the  tissue-level  mechanical  properties  that  we  measured  showed  a  small 
coefficient  of  variation  for  both  sexes,  approximately  20%  of  this  variation  was  explained  by  the 
size  of  the  tibia.  Female  tibiae  extended  the  correlations  that  were  observed  previously  for  males 

[12] ,  and  indicated  that  individuals  with  smaller  cross-sectional  size  (i.e.,  more  slender  bones) 
showed  a  larger  tissue-level  stiffness.  This  relationship  may  be  indicative  of  an  inherent  adaptive 
response  of  bone  to  modify  tissue-level  mechanical  properties  in  order  to  compensate  for  the 
smaller  size  [44].  This  coupling  suggested  that  osteoblasts  and  osteoclasts  are  capable  of 
coordinately  adapting  bone  morphology  and  tissue-level  mechanical  properties  so  that  the 
combination  of  these  traits  satisfies  mechanical  demands.  A  more  dramatic  coordinated 
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relationship  among  morphological  traits,  tissue-quality,  and  mechanical  function  has  been 
observed  previously  across  a  large  range  of  bones  from  many  species  [45]  and  for  inbred  mice 
[46-48],  bats  [49],  gulls  [50],  and  polar  bears  [51].  Although  this  coupling  appears  advantageous 
for  ensuring  that  an  adequate  whole  bone  stiffness  is  achieved  for  day-to-day  activity,  the 
disadvantage  is  that  the  tissue-quality  factors  that  tend  to  make  bone  stiff  also  tend  to  make  bone 
less  ductile,  less  tough,  and  more  damageable  [52],  These  latter  tissue-quality  factors  are 
important  during  extreme  loading  conditions  [53]. 

The  relationship  between  bone  size  (section  modulus)  and  body  size  (body  weight  x  tibia 
length),  which  was  established  previously  by  Ruff  [13]  and  was  based  on  the  resistance  of  bone 
to  bending  loads,  was  significantly  different  between  males  and  females  (Table  3,  Fig.  3). 
Similar  differences  were  observed  by  others  [6,  54-56].  Given  that  adaptive  factors  appear  to 
match  bone  size  to  body  size  during  growth  [39-42],  this  difference  suggests  that  these  adaptive 
factors  no  longer  work  the  same  way  when  men  and  women  reach  adulthood.  This  latent 
difference  in  bone  adaptation  may  be  due  to  differences  in  post-pubescent  endocrine  factors  that 
regulate  the  relative  amounts  of  bone  added  to  the  periosteum  versus  the  endosteum  [6].  Thus, 
adult  females  would  be  expected  to  accumulate  more  damage  under  intense  loading  compared  to 
males  and  this  may  be  a  contributing  factor  to  the  larger  stress  fracture  incidence  for  female 
athletes  and  military  recruits  [57,  58]. 

There  are  several  limitations  worth  addressing.  First,  the  difficulty  in  acquiring  tibiae  from 
healthy  men  and  women  in  the  20-40  year  age  range  limited  the  number  of  samples  for  this 
study.  Nevertheless,  this  did  not  affect  the  comparison  of  mechanical  and  morphological 
properties  between  sexes,  but  it  may  have  resulted  in  not  seeing  more  significant  correlations 
between  morphology  and  tissue-quality  (Table  4).  Second,  we  presented  only  measures  of 
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mechanical  properties,  and  have  not  explained  how  variability  in  these  mechanical  properties 
arises  from  variability  in  matrix  composition  and/or  micro-architecture.  Third,  tissue- 
damageability  was  based  on  measures  of  degradation  of  stiffness  and  relaxation.  Future  work 
will  determine  if  the  type,  distribution,  and  quantity  of  damage  is  different  for  females  and 
males.  In  addition,  our  measure  of  tissue-damageability  relates  to  the  progressive  accumulation 
of  damage  prior  to  failure  [19,  20],  and  may  not  reflect  how  damage  accumulates  under  long¬ 
term  fatigue  loading.  Fourth,  the  relationship  between  the  cross-sectional  size  of  the  tibia  and 
body  size  did  not  take  into  account  load  sharing  between  the  tibia  and  the  fibula.  Fifth,  our 
analysis  does  not  take  into  consideration  differences  in  the  amount  or  type  of  loading,  which 
have  been  shown  to  play  important  roles  in  the  development  of  bone  size  and  shape  [59].  Finally, 
our  analysis,  which  assumes  that  randomly  acquired  tibial  samples  provide  a  representative  view 
into  bone  biology,  does  not  take  into  account  lean  muscle  mass,  which  has  been  shown  to  be 
better  correlated  with  bone  size  and  shape  than  overall  body  weight  [60].  However,  these 
limitations  should  not  affect  the  similarity  in  mechanical  properties  between  males  and  females, 
or  the  correlation  between  bone  size  and  tissue-quality. 

The  results  of  this  study  provided  additional  insight  into  why  female  military  recruits  show  a 
greater  incidence  of  stress  fractures  compared  to  male  military  recruits  [57,  58].  The  undersized 
morphology  of  female  tibiae  combined  with  having  similar  tissue-level  mechanical  properties  as 
males  suggested  that  female  tibiae  are  overloaded  compared  to  males.  This  overloading  during 
intensive  training  would  be  expected  to  lead  to  increased  tissue-stresses  and  strains,  and 
subsequently  increased  in  situ  damage  accumulation  and  greater  risk  of  developing  a  stress 
fracture  [61-63].  Taken  together,  these  results  suggested  that  one  factor  contributing  to  the 
discrepancy  in  stress  fracture  risk  between  males  and  females  may  be  differences  in  the  way 
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bone  size  is  adapted  to  body  size  after  longitudinal  growth  has  ceased.  Thus,  fracture  risk 
reduction  would  benefit  from  having  a  better  understanding  of  the  factors  that  promote  or  inhibit 
bone  adaptation  [64-66]. 
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Table  1.  Monotonic  mechanical  properties  for  young-adult  females  and  males.  Data  were 
age-adjusted  and  shown  as  mean  +  standard  deviation. 


MECHANICAL  PROPERTY 

EEMALE 

MALE 

p-value 

Modulus,  GPa 

17.5  ±1.8 

17.1  ±  1.7 

0.5 

Yield  Strain 

0.008  ±0.0004 

0.008  ±0.0003 

0.7 

Yield  Stress,  MPa 

106.2  ±  10.0 

104.4  ±7.8 

0.6 

Post-yield  Strain 

0.023  ±0.006 

0.025  ±0.005 

0.3 

Eailure  Strain 

0.031  ±0.006 

0.033  ±0.005 

0.3 

Strength,  MPa 

133.8  ±8.4 

130.1  ±4.2 

0.2 

Energy-to-Eailure,  MPa 

3.1  ±0.7 

3.2  ±0.7 

0.3 
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Table  2.  Comparison  of  total  damage  accumulation  parameters  for  female  and  male 
cortical  bone.  Data  were  age-adjusted  and  shown  as  mean  +  standard  deviation. 


DAMAGE  PARAMETER 

EEMALE 

MALE 

p-value 

Stiffness  Degradation 

0.18  ±0.04 

0.17  ±0.03 

0.6 

Relaxation  Degradation 

1.89  ±0.08 

1.91  ±0.12 

0.6 
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Table  3.  Diaphyseal  cross-sectional  morphology  for  female  and  male  tibiae.  Data  were 
calculated  over  30-70%  of  the  total  tibial  length  and  shown  as  mean  +  standard  deviation. 


MORPHOLOGICAL  TRAIT 

FEMALE 

MALE 

p-value 

TIBIA  LENGTH,  mm 

33.9  ±2.8 

38.1  ±  1.9 

0.0001 

WIDTH AP,  mm 

26.5  +  2.7 

3L2±2.5 

0.0001 

WIDTHml,  nun 

20.8  +  2.4 

24.3  ±2.3 

0.0001 

CtAr,  mm^ 

249  +  40 

356  ±55 

0.0001 

J,  mm 

24776  +  7917 

51640 ± 15886 

0.0001 

J  /  rAp,  m^ 

1836  ±473 

3279  ±819 

0.0001 

J  /  Pml, 

2352  ±613 

4188  ±907 

0.0001 

J  /  L,  m^ 

722 ± 194 

1352  ±405 

0.0001 

J  /  A,  m^ 

97.1  ±  18.5 

142.7  ±24.9 

0.0001 

Sap,  1/mm  /kg 

14.1  ±4.1 

9.9  ±2.0 

0.002 

Sml,  1/mm^/kg 

11.1  ±3.6 

7.7  ±  1.5 

0.004 
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Table  4,  Correlation  matrix  between  morphological  traits  and  tissue-level  mechanical 
properties.  Combined  male  and  female  data  sets.  First  row  =  Pearson  eorrelation  eoeffieient. 


Seeond  row  =  p-value.  Damage  =  Stiffness  degradation. 


MODULUS 

STRENGTH 

PY  STRAIN 

ENERGY 

DAMAGE 

-0.13 

-0.15 

0.35 

0.34 

-0.12 

WIDTH Ap 

(0.5) 

(0.4) 

(0.05) 

(0.06) 

(0.5) 

WTFiTl^.  .X 

-0.13 

-0.17 

0.29 

0.28 

-0.03 

W IJJ 1  rlML 

(0.5) 

(0.4) 

(0.12) 

(0.13) 

(0.8) 

Pf  Ar 

-0.25 

-0.18 

0.31 

0.30 

-0.24 

(0.2) 

(0.3) 

(0.09) 

(0.1) 

(0.2) 

T 

-0.33 

-0.25 

0.27 

0.24 

-0.24 

J 

(0.07) 

(0.2) 

(0.14) 

(0.2) 

(0.2) 

T  /  T 

-0.35 

-0.26 

0.26 

0.22 

-0.26 

J  /  J-/ 

(0.05) 

(0.2) 

(0.2) 

(0.2) 

(0.2) 

T  /  Width  AT. 

-0.36 

-0.25 

0.22 

0.18 

-0.27 

J  /  VV  iU-LilAP 

(0.04) 

(0.2) 

(0.2) 

(0.3) 

(0.1) 

T  /  WidthATT 

-0.36 

-0.26 

0.24 

0.20 

-0.30 

J  /  VV  iU-LilML 

(0.05) 

(0.2) 

(0.2) 

(0.3) 

(0.1) 

c  . 

0.25 

-0.05 

-0.10 

-0.09 

0.43 

^AP 

(0.2) 

(0.8) 

(0.6) 

(0.7) 

(0.02) 

Qx  XX 

0.27 

-0.04 

-0.13 

-0.11 

0.44 

^ML 

(0.1) 

(0.8) 

(0.56) 

(0.6) 

(0.01) 
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Figure  1.  Females  show  an  age-related  change  in  tissue-level  properties  that  is  similar  to 
males,  (A)  modulus,  (B)  strength,  (C)  energy-to-failure,  and  (D)  damageability 
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Figure  2,  Changes  in  stiffness  and  relaxation  for  males  and  females.  A  plot  of  (A)  stiffness 
degradation,  Dstiff  and  (B)  relaxation  degradation,  Drelax,  for  each  diagnostic  cycle  in  the 
damage-accumulation  protocol.  Females  and  males  show  similar  curves.  Data  were  not  age- 
corrected. 
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Figure  3,  Variation  in  the  section  modulus  (J  /  WidthAp)  versus  body  size  (body  weight  x 
tibial  length)  for  adult  female  and  male  tibiae. 
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Figure  4,  Variation  in  mechanical  properties  as  a  function  of  diaphyseal  cross-sectional 
morphology.  Solid  lines  represent  the  regression  for  the  eombined  female  and  male  datasets. 
Dashed  lines  represent  the  regressions  for  males  and  females  separately. 
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INTRODUCTION:  Currently,  having  a  narrow  tibia  relative  to  body 
mass  is  one  of  the  best  predictors  of  stress  fracture  risk  and  fragility 
fractures  [1].  However,  the  reasons  are  not  fully  understood. 

Studies  of  bone  morphology  and  bone  quality  in  genetically  distinct 
inbred  mouse  strains  revealed  mice  with  slender  bones  had  increased 
mineral  content  compared  to  mice  with  more  robust  bones  suggesting 
that  bone  morphology  and  quality  might  be  biologically  coupled  to 
satisfy  mechanical  demands  imposed  by  weight  bearing  [2].  A  more 
strongly  coordinated  relationship  among  morphological  traits,  tissue- 
quality,  and  mechanical  function  has  also  been  observed  previously 
across  a  large  range  of  bones  from  many  species  including  bats,  gulls, 
and  polar  bears  [3].  Although  increased  mineral  content  may 
compensate  for  smaller  morphology  by  increasing  tissue  stiffness  and 
strength,  this  mineral  has  the  adverse  effect  of  increased  bone  brittleness 
and  tissue  damageability  under  fatigue  loading.  We  postulate  that  a 
similar  relationship  between  bone  morphology  and  quality  previously 
demonstrated  in  a  mouse  model  also  exists  in  the  human  skeleton.  To 
test  this  hypothesis,  we  quantified  cross-sectional  morphology,  tissue- 
level  mechanical  properties,  tissue  microstructure,  and  tissue 
composition  from  the  tibiae  of  young,  adult  males  and  females.  We  also 
tested  for  correlations  among  these  traits. 

METHODS:  Cross-sectional  morphology  and  bone  slenderness  of 
tibiae  from  14  female  donors  (age  22-46  yrs)  and  17  male  donors  (age 
17-46  yrs)  were  obtained  [4,5].  Mechanical  properties  were  measured 
from  8  cortical  bone  samples  (2.5mm  x  5mm  x  55mm)  machined  from 
the  diaphysis  of  each  bone  and  split  into  2  tests  groups  [4,5].  First, 
monotonic  failure  properties  were  assessed  by  loading  to  failure  in  4- 
point  bending.  Mechanical  properties  measured  were  modulus  (E), 
strength,  work,  and  post-yield  strain  (PYs)  as  a  measure  of  brittleness. 
Second,  tissue  damageability  was  assessed  using  a  fifteen-cycle  damage 
accumulation  protocol  in  4-point  bending  as  described  previously  [4,5]. 

The  density,  ash  content,  and  water  content  were  determined  for  each 
sample  retrieved  from  the  monotonic  tests.  Specimen  volume  was 
determined  using  Archimedes’  principal.  Submerged  weight,  hydrated 
weight,  dry  weight,  and  ash  weights  were  measured.  To  test  for  variation 
in  matrix  organization  in  males  and  females,  bone  microstructure  was 
assessed  for  each  sample  retrieved  from  the  damageability  tests.  For 
each  specimen,  digital  images  of  three  transverse  sections,  lOOpm  in 
thickness,  were  taken  at  lOX,  stitched  together,  and  traced  using  a  Tablet 
monitor  (WACOM).  Parameters  measured  include  porosity  and  the  area 
fractions  of  osteonal,  interstitial  (remodeled),  and  circumferential 
lamellar  (unremodeled)  tissues.  Both  vascular  canals  and  resorption 
spaces  were  counted  as  pores.  Osteonal  tissue  was  defined  as  a  lamellar 
region  with  a  Haversian  canal  completely  surrounded  by  a  cement  line. 
Data  from  individual  test  samples  were  averaged  for  each  donor. 

To  determine  if  males  and  females  show  a  similar  relationship 
between  tibia  cross-sectional  morphology,  bone  microstructure,  and 
tissue-level  mechanical  properties,  linear  regressions  were  performed 
while  taking  age  into  consideration  and  slopes  and  intercepts  between 
male  and  female  regressions  were  compared  (ANCOVA). 

RESULTS:  Since  no  differences  were  found  between  male  and  female 
regressions,  combined  male  and  female  datasets  were  used  to  compare 
tissue-level  mechanical  properties,  cross-sectional  morphology, 
composition,  and  bone  microstructure.  The  only  significant  correlations 
between  compositional  or  architectural  traits  and  tissue-level  mechanical 
properties  were  observed  between  ash  content  and  tissue  modulus  (p  < 
0.003),  strength  (p  <  0.02),  and  post-yield  strain  (p  <  0.01;  Fig.  lA). 
Measures  of  bone  size  (total  cross-sectional  area,  cortical  area,  ML  and 
AP  width,  cross-sectional  polar  moment  of  inertia)  were  negatively 
correlated  with  ash  content  (p  <  0.02).  Bone  slenderness  (polar  moment 
of  inertia  normalized  for  width  and  polar  moment  of  inertia  normalized 
for  tibia  length)  was  also  negatively  correlated  with  ash  content  (p  < 
0.02;  Fig.  IB).  Bone  size  and  slenderness  were  positively  correlated  with 
the  area  fraction  of  remodeled  tissue  (p  <  0.05).  Thus,  wider  (more 


robust)  tibiae  had  lower  mineral  content  and  increased  amount  of 
remodeled  tissue  and  narrow  (more  slender)  tibiae  had  increased  mineral 
content  and  decreased  amount  of  remodeled  tissue. 


Figure  1:  Correlation  between  ash  content  and  (A)  post-yield  strain, 
PYs  and  (B)  polar  moment  of  inertia,  J,  normalized  for  tibia  length,  L. 
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DISCUSSION:  The  current  study  confirmed  that,  similar  to  inbred 
mice,  the  tibia  diaphyses  of  young,  adult  females  and  males  show  a 
reciprocal  relationship  between  bone  morphology  and  tissue  quality,  as 
hypothesized.  Tissue  from  more  slender  bones  (narrow  relative  to  tibia 
length)  was  composed  of  material  having  higher  mineral  content  and, 
consequently,  increased  tissue  stiffness  and  strength,  but  reduced 
ductility.  Unlike  the  mouse  model,  this  relationship  was  subtle  and  only 
observed  when  combining  male  and  female  datasets  to  expand  the  range 
of  data.  The  current  study  revealed  a  general  biological  phenomenon  that 
has  only  been  observed  previously  by  looking  across  a  large  range  of 
bones  from  many  species.  Here,  bone  from  the  same  species  with  similar 
function  showed  that  more  slender  morphology  might  be  compensated 
by  increased  ash  content.  Although  this  coupling  appears  advantageous 
for  ensuring  that  an  adequate  whole  bone  stiffness  is  achieved  for  day- 
to-day  activity,  the  disadvantage  is  that  the  tissue-quality  factors  that 
tend  to  make  bone  stiff  also  tend  to  make  bone  less  ductile,  less  tough, 
and  more  damageable. 

This  is  an  intriguing  finding  because,  despite  dimorphic  growth 
patterns,  males  and  females  have  similar  tissue-level  mechanical 
properties.  Thus,  both  sexes  may  share  a  common  underlying  biological 
(adaptive)  mechanism.  This  data  suggest  that  osteoblasts  and  osteoclasts 
are  capable  of  not  only  adapting  bone  morphology,  but  also  modulating 
tissue-level  mechanical  properties  so  that  the  combination  of  traits 
satisfies  mechanical  demands  [2,3].  Therefore,  the  phenomenon  that 
slender  bones  are  constructed  differently  with  material  level  variation 
that  ultimately  leads  to  more  damageable  material  than  larger  bones  may 
help  explain  why  individuals  with  more  slender  bones  are  at  increased 
risk  of  stress  fractures  early  in  life  and  fragility  fractures  later  in  life. 
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